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Chapter  1 

Introduction  and  Summary 


A  sensitive  superconducting  six-axis  accelerometer  (SSA)  has  been  developed 
for  application  in  gravity  survey  and  inertial  guidance.  A  breadboard  SS.4  has  been 
constructed  and  operated  to  verify  the  design  concepts  of  the  device.  The  construction 
of  the  breadboard  SSA  followed  closely  the  original  design  proposed  by  Paik  [d,ll].  Based 
on  the  experience  obtained  with  the  breadboard,  an  improved  model  of  the  SSA  has  been 
designed  and  is  under  construction. 

The  design  involves  a  single  levitated  superconducting  proof  mass  of  a  special 
geometry  whose  motion  is  monitored  in  six  degrees  of  freedom.  The  proof  mass  has  the 
shape  of  c.n  inverted  cube,  three  orthogonal  square  plates  bisecting  each  other,  and  is 
made  of  niobium  (Nb).  It  is  levitated  on  magnetic  fields  produced  by  persistent  currents 
in  24  superconducting  levitation /feed back  coils.  Its  motion  in  three  linear  and  three 
angular  degrees  of  freedom  is  sensed  by  24  superconducting  sensing  coils,  which  form  six 
independent  inductance  bridges.  The  six  bridges  are  driven  at  six  different  frequencies 
and  are  detected  by  a  single  SQUID,  whose  output  is  demodulated  to  recover  the  six 
acceleration  signals.  The  persistent  currents  in  the  levitation/feedback  coils  are  adjusted 
to  null  the  six  inductance  bridges  and  the  SQUID  demodulated  signals  are  fed  back  to 
the  corresponding  levitation/feedback  coils  to  keep  the  proof  mass  at  the  null  position. 

The  coil  forms  have  been  machined  out  of  a  titanium  alloy,  TiV4Al6,  and  all 
the  coils  were  v.-ound  initially  with  niobium  wire.  Niobium  was  chosen  because  of  its 
high  value  of  //d  so  that  the  coils  can  be  operated  in  the  Type-I  region  where  currents 
are  more  stable.  All  the  parts  fit  inside  a  10.16  cm  titanium  cube. 

Two  problems  were  encountered  with  this  design.  First,  the  titanium  alloy 
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used  for  the  coil  forms  had  an  unusually  low  superconducting  transition  temperature, 
4.5  K.  Since  the  magnetic  fields  should  be  able  to  return  through  the  coil  forms,  we  had 
to  operate  the  SSA  at  an  elevated  temperature  near  5  K  to  keep  the  coil  forms  from 
becoming  supeicOiiducting.  In  the  new  model  of  the  SSA,  a  different  alloy  of  titanium  is 
used  whose  transition  temperature  is  below  1.5  K.  Another  problem  that  has  plagued  i  lie 
progress  of  our  research  is  the  brittleness  of  the  thin  niobium  wires  at  low  temperature'-. 
After  a  long  battle  to  overcome  the  problem  of  wire  breaking  upon  thermal  cycling.  \v< 
resorted  to  a  compromise  solution  of  winding  the  coils  with  niobium-titanium  (.Nb  l  ii 
wire.  Although  mechanically  strong,  the  Nb-Ti  wire  has  to  be  operated  in  a  Type-Il 
region,  thereby  increasing  the  drifts  in  the  superconducting  circuits.  In  the  improved 
model,  we  have  gone  back  to  Nb  wire  to  eliminate  this  problem. 

The  breadboard  SSA  has  been  operated  completely.  Levitation  and  balanre* 
have  been  achieved  in  all  six  degrees  of  freedom.  The  outi)ut  of  a  single  SQl’lD  has  b('en 
successfully  demodulated  to  recover  the  six  acceleration  signals.  The  external  rf  noise, 
which  had  been  a  main  concern  for  operating  the  bridge  s^asing  circuits  at  the  input  of 
a  SQUID,  has  been  filtered  sufficiently  by  means  of  rf  shielded  transformers. 

In  addition  to  constructing  an  improved  model  of  SS.‘\.  further  tests  will  be 
carried  out  on  the  breadboard  SSA  under  a  new  .\ir  Force  crnitract.  In  particular, 
the  effect  of  oscillator  noise  and  the  source  of  any  excess  low  frequency  noise  must  be 
investigated.  The  feedback  circuits  need  to  be  developed  to  improve  the  dynamic  range 
of  the  device. 

The  theoretical  sensitivity  of  the  SS.'\  is  10~'^g  for  the  linear  and 

10“®  arcsec  sec“^  for  the  angular  acceleration.  The  SS.-\  of  the  new  design  will 

be  integrated  with  the  three-ajds  superconducting  gravity  gradiomcter  (SCG)  being  de¬ 
veloped  under  a  NASA  contract  It  is  pnvision^d  that  the  combined  SGG/SSA  package 
will  be  flown  in  the  Earth  orbit  for  a  gravity  mapping  mission  in  the  late  1900's.  The 
superconducting  gravity  gradiometer/acceleromeier  system  with  such  high  sensitivity, 
if  successfully  developed,  will  find  many  important  applications  in  gravity  survey  and 
inertial  guidance. 
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Chapter  2 


Theory  of  a  Superconducting 
Six- Axis  Accelerometer 

2.1  Introduction 

A  "six  axis"  ar roloromelpr  measures  the  linear  and  angular  arc rderal ion  in  all 
six  degrees  of  freedom  at  the  same  point  in  spare*  lime*.  A  superrondurting  six-axis 
arreleroreeter  opc'rates  at  lic|iiid  helium  temperatures  and  is  able  to  take  advantage  of 
the  low  thermal  noise  and  the  extremely  high  stability  of  materials  and  currents,  whirh 
is  a  direct  conseciuence  of  cryogenic  temperatures.  The  six-axis  accelerometer  uses  the 
properties  of  quantized  magnetic  flux  to  magnetically  levitate  a  single  superconducting 
proof  mass  with  extreme  stability.  Additional  superconducting  circuits  monitor  and 
control  the  position  cjf  the  proof  m.iss.  High  sensitivity  is  achieved  through  the  use*  of  a 
.SQUID  (Superconducting  QUanturn  Interference  Device)  to  amplify  the  displacements 
of  the  proof  mass  away  from  equilibrium. 

Like  the  accelerometers  in  the  gravity  gradiometer.  the  superconducting  proof 
mass  in  the  six-axis  accelerometer  (SSA)  is  free  to  respond  to  accolcraticn.  However,  the 
proof  mass  in  this  case  is  free  to  respond  to  accelerations  in  all  six  degrees  of  freedom.  In 
the  accelerometers  in  the  gravity  gradiometer,  the  proof  ma.ss  position  is  confined  by  a 
mechanical  spring  and  the  accelerometer  responds  to  arc('lcratlon  in  iiuerse  proportion 
to  the  stiffness  of  this  mechanical  spring.  In  the  six  axis  .accelerometer,  the  proof  mass 
is  surrounded  by  quantized  magnetic  flux  whirh  confines  the  motion.  Each  degree  of 


tVecdoiii  resj)oiui,s  to  arrolt’ralioii  in  invi'isc  [)roporlion  to  tio-  >lifriK'ss  of  a  ■■niafiin  tir 
spring”. 


Those  “niagiietic  springs”  ar*-  goiiorat<'ii  by  the  sup^'roond net ing  roilv  placet! 
in  close  proximity  to  the  proof  mtiss.  I’ho  niobiii:ii  proof  mass  e\(  ludo.•^  the  magnet k 
flux  from  these  coils  due  to  the  Meissner  effect  and  any  motion  of  the  proof  rlla'^  for'  ■ 
flux  to  redistribute  itself  within  the  superconducting  (ireuitrv.  This  transfer  of  en.  i'j  v 
produces  a  restoring  forc<'  which  opposes  this  change  and  creates  the  electrical  etpiival.  at 
to  a  mechanical  spring. 

Each  degre(' of  fre<'dom  in  iheSS.A  isconfinc'ii  and  (oiitrolled  by  ;in  indeperident 
superconducting  circuit  made  ii[)  of  four  levitation  coils.  1  lie  stiffness  in  each  degiee  of 
freedom  is  controlled  by  I  hese  circuits.  I'he  t hree  <  irenits  which  control  the  linear  degrees 
of  freedom,  jrrovide  the  le’.itation  forcr-  necessary  to  susi)end  the  proof  mass  against  the 
pull  of  gravity.  'I'he  thrr'e  renutining  circuits  proside  feedback  control  for  the  ang'ila. 
degrees  of  freedom.  In  all.  twenly  foiir  levitation  inductaiKes  surround  the  proof  ma.ss. 

•An  additional  twenty-four  srnising  inductances  suriound  the  [iroof  ma.ss.  aitd 
are  used  to  detect  ttie  motion  of  the  proof  imuss.  W'Ik'ii  an  ideration  is  ajrplitnl  to  the 
.S.SA.  the  proof  mass  is  disj)lace<l  in  in\a>rse  proportion  to  the  '•iitiness  of  the  "magncaic'' 
spritigs  wliich  surround  th<'  proof  ma.ss.  I)<'t<>ct if)n  of  this  disiibuement  is  [xissible  since 
tlte  .Mei.-.sn(’r  <'ffect  forces  the  inductance  of  the  smising  coils  to  chang''  in.  proportion 
to  the  displacenuMit  of  tin*  jrroof  mas.s.  Each  degrei'  of  freedr)m  i.s  imniitored  by  four 
sensing  coils  connected  togethi'r  to  form  a  superconducting  bridge  (irruit.  I  he  sensing 
coils  which  niiike  u[)  each  trridge  circuit  are  properly  selected  so  that  the  output  of  each 
liridge  is  de|)endent  on  onlv  fine  degree  of  freefbun.  I'he  outputs  of  all  tiu'  bridge  ciia  uit^ 
.are  connected  in  serif's  to  a  SQEII)  amplifif'r.  which  f.uiverts  tho  signal  to  an  output 
voltage.  Eac  h  bridge  is  driven  by  ,tn  indepf'ndent  oscillating  (  urrent.  .and  a  misbalancf’  of 
any  one  of  the  briilge  ci nulls  will  ajijiear  at  the  out  [nit  tif  tin'  E ID  ami>lifier.  Rf't au  f'ry 
of  the  original  disiihu  I’lneni  of  the  proof  ma.ss  in  f'ach  df'gree  of  freedom  is  ai  cotnplisl,.  d 
through  the  use  of  six  lock  in  amplifn'rs.  each  iruiependent ly  tuned  to  tin'  correct  carrier 
fn'fpiency. 

.A  sirnjilifif'd  Vf'rsifin  of  Inis  detection  tedinifpie  n  -.hown  in  Eig.  'd.l.  for  tu. 
a(  celeromef ef  w  itli  one  df'gref'  of  freedom.  .As  tlif'  prool  lu.i.^^  is  displaced  from  ('(juilih 
nuni  in  the  [msilive  r  direction,  the  sensing  induit.mces  /_  and  /  j  dec n-ase  and  the 
sf'nsing  I  ml  u  r  t  aiK  cs  /.  |  .mil  /,  i  increase.  This  change  m  ■.eiising  induit.UKf'  unbalances 


S  upc  rcond  acting 

l<.,n.Uc  Sprin, 


Bridg* 

Circuitry 


SQUID 


o(iu)  - -  x{w) - —  L(w)  t(u/)costi4t  - -  l^iu) 

Acceleration  Displacement  inductance 

Modulation 


F'iguro  2.1.  .A  schomatir  illustrating  tho  principle  of  a  s)i[)crcoiHlu(  ting 
bridge  accoloratioa  transducer. 


the  bridge  circuit  and  an  oscillating  current  proiiortional  in  amplitude  to  the  misbalance 
appears  across  the  output  coil.  This  current  is  then  amplified  and  converted  to  an  out¬ 
put  voltage  through  the  use  of  a  SQUID  amplifier.  Demodulation  through  the  use  of  a 
lock-in  amplifier  then  recovers  the  original  misbalance  signal. 

In  this  section,  we  have  pre.sented  a  very  brief  di'scription  of  the  operation  of 
the  SSA.  In  the  rest  of  this  chapter,  we  will  expand  upon  this  descrijition.  We  will  first 
present  a  general  description  of  the  dynamicaJ  equations  of  the  accelerometer.  This  will 
be  follow'ed  by  a  description  and  analysis  of  the  levitation  and  simsing  circuitry.  We  will 
then  combine  this  analysis  and  derive  the  transfer  functions  for  the  SS.A  which  relate 
the  external  acceleration  of  the  housing  to  the  output.  In  the  final  section,  we  disciis> 
the  potential  sensitivity  of  the  SSA.  Throimhout  this  chapter,  unless  otherwise  statod. 
we  will  employ  sumr.iation  notation,  wlmre  any  two  like  indices  will  denote  a  sum  over 
each  of  the  three  axes,  t.  y,  and  c. 
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2.2  Dynamics  of  a  Superconducting  Six- Axis 
A  ccelerometer 

An  accelerometer  must  measure  the  acceleration  of  itself  relative  to  an  inertial 
reference  frame.  In  a  six-axis  accelerometer,  it  is  necessary  to  simultaneously  measure 
the  three  linear  and  the  three  angular  accelerations  applied  to  the  accelerometer.  In 
actuality,  the  position  and  orientation  of  the  proof  mass  is  measured  relative  to  t;.  ' 
accelerometer  housing,  and  the  acceleration  of  the  proof  mass  housing  is  inferred  from 
the  dynamics  of  the  proof  mass.  Thus,  it  is  necessary  to  clearly  understand  how  th^" 
dynamics  of  the  accelerometer  ran  he  used  to  derive  the  applied  acceleration  from  a 
knowledge  of  the  position  and  orientation  of  the  proof  mass. 

In  an  accelerometer  with  only  one  degree  of  freedom,  the  equation  of  motion 
that  relates  the  acceleration  of  the  accelerometer  housing  to  the  position  of  the  proof 
mJLSS  is  relatively  simple.  In  the  SSA,  it  is  necessary  to  contend  with  both  the  linear  and 
angular  degrees  of  freedom.  This  additional  angular  freedom  complicates  the  description 
greatly  and  we  wiU  first  consider  the  dynamics  of  a  single-axis  accelerometer. 

We  will  then  derive  the  dynamical  equations  for  a  rigid  body  in  rotating,  ac¬ 
celerated  reference  frame,  i.e.,  the  dynamical  equations  of  the  proof  mass  relative  to  the 
accelerometer  housing  which  is  rotating  and  accelerated.  These  dynamical  equations 
will  consist  of  three  linear  equations  of  motion  and  three  Euler  equations.  The  three 
linear  equations  will  directly  relate  the  measurement  of  the  position  of  the  proof  mass 
to  the  acceleration  of  the  accelerometer  housing.  The  three  angular  equations  will  be 
more  complicated.  Fiuler’s  equations  provide  a  description  of  the  angular  acceleration  of 
a  body.  It  will  be  necessary  to  parametrize  the  orientation  of  the  proof  mass  in  terms 
of  a  particular  set  of  angles,  and  combine  the  kinematic  ecpiations  with  the  dynamical 
equations  to  derive  a  set  of  angular  equations  of  motion  that  relate  the  orientation  of 
the  proof  mass  to  the  angular  accelerations  of  the  proof  ma.ss  housing. 

2.2.1  Dynamics  of  a  Single-Axis  Accelerometer  in  a  Rotating,  Accel¬ 
erated  Reference  Frame 

In  this  section,  we  will  be  concerned  with  two  referemee  franu's:  the  referenct' 
frame  attached  to  the  accelerometer  housing  from  which  the  proof  ma.ss  positiim  is 
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Figure  2.2.  The  single-axis  accelerometer  proof  mass  coordinate  frame. 

measured,  and  an  inertial  reference  frame  in  which  the  accelerometer  housing  is  rotating 
and  accelerating.  In  order  to  derive  Lagrange’s  equations,  we  will  work  in  component 
notation.  This  is  not  as  elegant  as  working  in  vector  notation,  but  it  is  powerful.  After 
we  have  derived  Lagrange’s  equations,  we  can  return  to  the  elegance  of  vector  notation. 

The  proof  ma.ss  of  a  single-axis  accelerometer  is  shown  schematically  in  Fig. 
2.2.  The  unprimed  coordinates  are  taken  as  fixed  to  the  accelerometer  housing,  and  "he 
primed  coordinates  are  fixed  to  an  inertial  reference  frame.  The  coordinates  of  the  cen  .  r 
of  mass  of  the  proof  mass  in  the  inertial  reference  frame  are  defined  by 

•^i  =  (2.1) 

whore  is  the  rotation  matrix  descril)ing  the  orientation  of  the  accelerometer  housing 
with  respect  to  the  inertial  reference  frame,  is  the  position  of  the  accelerometer 
housing  in  the  inertial  reference  frame,  and  is  the  position  of  the  proof  mass  body 
relative  to  the  accelerometer  housing. 

In  order  to  formulate  the  dynamics  of  the  accelerometer,  we  must  compute  the 
velocity  of  the  proof  mass  in  the  inertial  reference  frame.  This  is  a  function  of  7?^^, 

,  and  .  Once  we  know  this  velocity,  we  can  write  down  the  Lagrangian  for  the 
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accelerometer  and  solve  for  the  equations  of  motion.  Differentiating  Eq.  (2.1)  with 
respect  to  time  gives  the  components  of  the  velocity  of  the  center  of  mass  in  the  inertial 
reference  frame, 


u'  =  itHr^  +  RUH’  +  r! 


H-P 


uP 


(2.2) 


This  can  be  written  in  the  standard  way  by  grouping  terms  .so  that  all  coristit  aeiit 
components  can  be  viewed  in  the  inertial  reference  frame.  We  have 


=  RliRUf'R 


H.P  ,  nh-P  ,  -H 
jl^l  +  ^ik^k  +  ^  • 


(2.3) 


The  second  term  in  this  equation  is  the  velocity  of  the  proof  mass  relative  to  the  housing 
as  viewed  in  the  inertial  reference  frame,  and  the  third  term  is  just  the  velocity  of  the 
housing.  The  first  term  is  the  product  of  two  terms:  the  first,  RIIR{'j  ' ,  is  related  to  the 
angular  velocity  of  the  housing,  and  the  second,  Rj^irf,  is  just  the  position  of  the  proof 
mass  relative  to  the  housing  as  seen  in  the  inertial  reference  frame.  From  Fq.  (2.3)  we  see 
that  RflR^j  *  maps  the  position  of  the  proof  mass  in  the  rotating  frame  into  its  velocity 
as  seen  from  the  inertial  reference  frame.  Furthermore,  R^^R^j  '  is  antisymmetric  and 
related  to  the  components  of  angular  velocity,  u).  In  the  inertial  reference  frame,  the 
components  of  o)  satisfy  the  equation  [5]: 

R'/^R^,"  (2.4) 

This  is  the  usual  w.  Substituting  Eq.  (2.4)  into  Eq.  (2.3).  we  have 

v'  =  ujj(,jkRi{r[‘  +  Rilrl  +  .  (2.-5) 


We  may  now  rewrite  this  in  vector  notation: 

v'  —  UP  X  . 


(2.6) 


where  is  the  velocity  of  the  proof  mass  relative  to  the  liousiiig,  and  is  the  velocity 
of  the  housing.  We  have  derived  this  standard  result  from  a  component  formalism.  I'his 
is  the  same  technique  we  will  use  to  derive  Lagrange's  equations  of  motion. 

The  Lagrangian  for  the  dynamics  of  an  accelerometer  in  a  rotating,  accelerated 
reference  frame  is 


(2.7) 
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where  V(r^)  is  the  potential  energy  internal  to  the  accelerometer,  and  is  the 

gravitational  potential  energy.  The  internal  potential  energy,  V,  includes  the  energy  of 
the  mechanical  spring  and  the  electromagnetic  energy  of  the  superconducting  circuitry. 

The  generalized  linear  momentum  from  this  Lagrangian,  as  measured  with  re¬ 
spect  to  the  accelerometer  housing,  is 


pi^=  ^  ^  ■ 

dn  94 

The  generalized  force,  as  measured  with  respect  to  the  experimental  housing,  is 
dC 


fk  = 


dr^ 


dv] 

—  m  -  v'l  ~  m 

drC  ' 


d<^B(x')  dx'  c)V'(r^) 


dx'  drl 


arf 


=  (^R^i  rf  +  rf  +  r")  -  tuRU 


d^E(x')  dV(r^ 


dx[ 


(2.8) 


(2.9) 


The  time  derivative  of  the  generalized  linear  momentum  can  now  be  computed; 


=  +  R^r  +  r'/) 

^  mR%{R>lrr^2RX+ li'Jf I" +  <=!')■ 


(2,10) 


Lagrange’s  equation  sets  the  generalized  force  equal  to  the  time  rate  of  change  of  the 
generalized  momentum: 

d^Ei^')  1  dVir^) 


R^.iR^.rf  +  2R^ff  +  fljji-r  +  r«)  -  -R^k 


fUr  X.  pW-  P  . 


dx',  ^  drt 


(2.11) 


Using  the  fact  that  is  orthogonal,  R^[  =  R^-  ',  and  rearranging  terms  gives 

1  OVir'') 

rl  - - =  af , 

dri' 

where  is  the  external  acceleration  applied  to  the  accelerometer: 

d^Ei^')  \ 


ai--R^r  (rj'+2«;;r;  +/?«r/'  + 


2P 


Hi.p  .  pW.P 


dx'. 


(2.12) 


(2.13) 
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The  external  acceleration  applied  to  the  accelerometer,  ,  can  be  written: 

d^Eix')  \ 

+  - ^  R^,rf  -  Ri'~  R^,rf.  (2.14) 

dx'^  J 

This  can  be  written  more  clearly  in  terms  of  the  angular  velocity  of  the  housing,  uj.  We 
have  already  seen  the  components  of  u>  in  the  inertia!  reference  frame,  Eq.  (2.4).  The 
components  of  u,  as  viewed  from  the  rotating  housing,  are  defined  as 


w,- 


(2.15) 


In  analogy  to  Eq.  (2.4),  we  find  that 

Rkf  R^t  =  (2.16) 

can  similarly  be  replaced.  By  differentiating  Eq.  (2.16)  with  respect  to  time 
and  applying  Eq.  (2.16)  to  the  result,  we  find 

Rkj  Rjt  —  kjm^mil  T  ^i^kih  (2.1  (  ) 


Substituting  Eqs.  (2.16)  and  (2.17)  into  Eq.  (2.14),  we  find 

d^E{x') 


af  =  -iZg-' 


Tj  + 


-  2u>jrf(kjl  -  ikjm^}(mtl^,rf  -  ikjm’^iT^' .  (2.18) 


dx'j 


If  we  rewrite  this  in  vectorial  notation,  the  result  is  clear: 


-  V$£  -  2u?xf^-u>x(u>x  r^)  -  u>  X  .  (2.19) 


The  first  term  of  the  right  hand  side  is  the  acceleration  of  the  housing,  the  second  term 
is  the  gravitational  acceleration,  the  third  term  is  the  Coriolis  acceleration,  the  fourth 
term  is  the  centrifugal  acceleration,  and  the  fifth  and  final  term  couples  the  angular 
acceleration  to  the  linear  acceleration  in  proportion  to  the  displacement  of  the  proof 
mass  away  from  the  center  of  the  accelerometer  housing. 


2.2.2  Dynamics  of  a  Rigid  Body  in  a  Rotating  Accelerated  Reference 
Frame 

In  the  SSA,  it  is  necessary  to  contend  with  both  the  linear  and  angular  degrees 
of  freedom.  This  additional  angular  freedom  complicates  the  description  greatly.  We  will 
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Figure  2.3.  Schematic  of  the  six-axis  accelerometer  proof  mass. 

be  dealing  with  three  specific  reference  frames:  an  inertial  reference  frame,  a  reference 
frame  attached  to  the  accelerometer  housing,  and  a  reference  frame  attached  to  the  proof 
mass.  Since  all  measurements  of  the  position  and  orientation  of  the  proof  mass  are  made 
with  respect  to  the  accelerometer  housing,  we  will  derive  all  of  our  equations  of  motion 
with  respect  to  that  coordinate  system. 

We  will  begin  the  derivation  of  the  dynamical  equations  by  describing  the  po¬ 
sition  and  orientation  of  the  proof  mass  in  terms  of  our  three  specific  reference  frames. 
After  this  is  done,  we  will  derive  the  linear  and  angular  velocity  of  the  proof  mass  in  iner¬ 
tial  space.  This  will  allow  us  to  construct  the  Lagrangian  for  SSA.  After  the  Lagrangian 
is  constructed,  we  will  derive  our  linear  and  angular  equations  of  motion. 

The  proof  mass  of  the  six-axis  accelerometer  is  shown  in  Fig.  2.3.  The  unprimed 
coordinates  are  taken  as  fixed  in  the  body  of  the  proof  mass,  the  primed  coordinates  arc 
fixed  to  the  accelerometer  housing,  and  the  double  primed  coordinates  are  fixed  to  an 
inertial  reference  frame. 

In  order  to  understand  the  dynamics  of  the  proof  mass,  we  begin  by  first  relating 
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the  components  of  a  single  point,  P,  as  seen  in  the  three  coordinate  systeins.  !n  the 
reference  frame  fixed  to  the  proof  mass,  the  coordinates  of  P  are,  by  definition. 


(2.20) 


This  same  point,  as  seen  in  the  primed  coordinate  system  attached  to  the  accelerometer 
housing,  has  coordinates 


I'i  =  + 


(2.21) 


where  is  the  rotation  matrix  which  describes  the  orientation  of  the  proof  mass  with 
respect  to  the  accelerometer  housing,  and  is  the  displacement  of  the  proof  mass  away 
from  the  center  of  the  housing.  The  coordinates  of  P,  as  seen  in  the  inertial  reference 
frame,  are 


R 


H 


H 


(2.22) 


where  is  the  rotation  matrix  which  describes  the  orientation  of  the  housing  of  the 
accelerometer  with  respect  to  the  inertial  reference  frame,  and  is  the  position  of  the 
housing  in  the  inertial  reference  frame. 

In  order  to  construct  the  Lagrangian  which  describes  the  dynamics  of  the  SSA, 
we  must  compute  the  angular  velocity  and  the  linear  velocity  of  the  proof  mass  relative 
to  inertial  space.  This  velocity  will  be  a  function  of  R'^  ,  and  .  Once  this  is 
done,  we  will  be  able  to  solve  for  the  equations  of  motion. 

We  are  free  to  choose  the  axes  of  the  coordinate  system  fixed  to  the  proof  mass 
to  lie  along  the  principal  body  axes  of  the  proof  mass,  and  free  to  choose  the  center  of 
mass  to  lie  at  the  origin  of  this  coordinate  system.  Thus  the  coordinates  of  the  center 
of  mass  in  the  inertial  reference  frame  are 


r"  =  R 


IKP 


f  T 


H 


(2.23) 


where  we  have  set  x  =  0  in  f’q.  (2.22).  This  is  exactly  the  same  equation  a.s  in  Section 
2.2.1,  wheie  we  described  the  dynamics  of  a  single-axis  accelerometer  in  a  rotating  ac¬ 
celerated  reference  frame.  Following  the  same  procedure  as  in  Eqs.  (2.1)  through  (2..5). 
the  velocity  of  the  center  of  mass  of  the  proof  mass  is 


=^j  +  Rik^k  + • 


(2.24) 
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where  is  the  angular  velocity  of  the  housing  in  inertial  space,  and  has  components 
which  satisfy  the  equation: 

(2.25) 


We  may  now  rewrite  Eq.  (2.24)  in  vector  notation: 


V  J.  -L 

V  —  (jJ  XT  -j-  V  i-tj  , 


(2.26) 


where  u  ”  is  the  velocity  of  the  center  of  mass  of  the  proof  mass  relative  to  the  housing, 
and  is  the  velocity  of  the  housing. 

Next,  v;e  must  compute  the  angular  velocity  of  the  proof  mass  with  respect  to 
inertial  space.  In  this  case  the  displacements  of  the  coordinate  origins,  and  do 
not  affect  the  transformation  of  a  vector.  In  fact,  if  /  is  an  arbitrary  vector  fixed  to  the 
accelerometer  proof  mass,  then  we  may  write 


/"  =  R^,R^,fk. 


(2.27) 


Difierentiating  this  vector  with  respect  to  time  must  give  us  ihe  angular  velocity.  We 


and  this  may  be  rewritten: 


(2.28) 


/"  =  R^R^  RimRLJk  +  Rki~ 


Substituting  Eq.  (2.27),  we  find 


/"  =  {M^Ry  +  R^'R^Xk'Rkr')fl'- 


(2.29) 


(2.30) 


Because  /  is  a  vector  fixed  to  the  accelerometer  proof  mass,  its  time  derivative  must 
obey  the  equation 

/:'  =  c.X7r-  (2.31) 

Identifying  the  two  coefficients  of  //'  in  Eqs.  (2.30)  and  (2.31),  we  have 

UjlU,’]  =  +  R^,R^rRrk'Rki~'-  (2.32) 

The  first  term  on  the  right  hand  side  of  this  equation  is,  from  Eq.  (2.4), 

R»R^,~'  =ufuji-  (2.33) 
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The  second  term  on  the  right  hand  side  of  Eq.  (2.32)  is  more  complicated.  The  two 
middle  factors,  can  be  w'ritten,  from  Eq.  (2.4): 

(2.34) 


Substituting  these  two  results  into  the  equation  for  the  angular  velocity  of  the  proof 
mass,  we  find 

"I"  ^ij^jrk^r  ^kl  '  (2.3  )) 

Recognizing  that  the  components  of  the  Levi-Civita  tensor  must  remain  the  same  in  all 
Cartesian  coordinate  systems,  this  can  be  written: 


Thus,  the  components  of  the  angular  velocity  of  the  proof  mass,  u>,  are 


(2..36) 


U)'J  =  +  ^]*k^'k  ■> 


(2.37) 


where  is  the  angular  velocity  of  the  accelerometer  housing,  and  R^ is  the  angular 
velocity  of  the  proof  mass  with  respect  to  the  housing,  as  seen  in  the  inertial  reference 
frame. 

We  can  now  compute  the  Lagrangian  for  the  six-axis  accelerometer: 

C  =  2^^^"  '  2^  •  I  ■  w  -  V(R^ ,r^)  -  in4>(x''),  (2.38) 

where  V{R^ is  the  potential  energy  internal  to  the  accelerometer,  (t>(x")  is  the 
gravitational  potential  energy  of  the  proof  mass,  m  is  the  mass  of  the  proof  mass,  and 
I  is  the  moment  of  inertia  tensor  for  the  proof  mass.  The  internal  potential  energy,  V^, 
includes  the  electromagnetic  energy  of  the  superconducting  circuitry. 


Linear  Dynamical  Equations 


We  first  compute  the  dynamical  equations  for  the  linear  degrees  of  freedom. 
The  generedized  linear  momentum  from  the  SSA  Lagrangian,  as  measured  with  respect 
to  the  accelerometer  housing,  is 


(2..39) 
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The  generalized  linear  force,  as  measured  with  respect  to  the  accelerometer  housing,  is 
/'■  -  - 

d^E(x"}  dx''  dV{R^,r^) 


dv" 

J  // 

=  m  -  V:  -  m 


dr^  ' 


dx'f  dr^ 


=  mk^,  [k^ivf  +  R^,rf  +  ff )  -  m/?«  ^  •  (2.40) 

£)x, 

This  analysis  is  identical  to  that  in  Section  2.2.1.  Setting  the  time  derivative  of  the 
generalized  momenta  equal  to  the  generalized  force,  and  following  the  same  analysis  as 
in  Eqs.  (2.3)  through  (2.8),  we  find 


1  OT(ii^r^) 

^r+- 


dri: 


(2.41) 


where  is  the  external  acceleration  applied  to  the  accelerometer: 


„E  _  r,H-‘  "H  , 

Q).  =  -Kkj  Xj  + 


d^E{x'') 


n  iH  -P  Itl  P  ■  IH  P 

-  2ulj"r,  (kjl  -  Ckjm^^  J  fmt/W  ,  r,  -  ffcjmW  ,  r;  , 


tip 


JP^P 


dx'l 


(2.42) 

where  are  the  components  of  the  angular  velocity  of  the  housing,  as 

viewed  from  the  rotating  housing.  If  we  rewrite  this  in  vectorial  notation,  the  result  is 
clear: 


-  V^E  -  2  01^^  X  X  X  j  -  X  r^.  (2.43) 

This  result  is  identical  to  Eq.  (2.19).  The  first  term  of  the  right  hand  side  is  the  accel¬ 
eration  of  the  housing,  the  second  term  is  the  gravitational  acceleration,  the  third  term 
is  the  Coriolis  acceleration,  the  fourth  term  is  the  centrifugal  acceleration,  the  fifth  and 
final  term  couples  the  angular  acceleration  to  the  linear  acceleration  in  proportion  to 
the  displacement  of  the  proof  mass  away  from  the  center  of  the  housing.  This  equation 
forms  the  cornerstone  for  an  iinderstanding  of  the  linear  acceleration  of  the  center  of 
mass  of  the  six-axis  accelerometer. 


Angular  Dynamical  Equations 

The  angular  equations  of  motion  can  now  be  computed.  In  order  to  correctly 
compute  the.se  equations,  we  must  incorporate  six  constraint  equations  that  insure  that 
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the  state  variables,  are  varied  in  such  a  way  that  the  matrix  remains  an  or¬ 
thogonal  transformation.  We  can  do  this  by  using  the  method  of  Lagrange  multipliers, 
whereby  we  add  our  constraint  equations  to  the  Lagrangian.  We  find 

^  •/  -  uj  -  -  m4>{x")  +  A;„,  -  6im)  ,  (2.41) 

where  is  symmetric.  This  will  insure  that  the  matrix  R^  remains  orthogonal,  i.e., 

«p-'5/m  =  0.  (2.45) 


The  generalized  angular  momentum  from  this  Lagrangian  can  now  be  computed: 

,P 


p'T  = 


dC 


OR 


bP 


UUJ  L  ijU)  I  /  rr  n  V 

t'l  .  _  dP  t'i  ( .  1  dH  . 

q  ~  ^kl  ~~~  ^“-7  u  ) 

d  R\ 


dR‘ 


du'^ 

dkf^ 


bP 

•j 

H  ,  ,P 
q  +‘^9 


). 


(2.46) 


where  we  have  made  the  substitution 


//-■ 

Iv 


Jtt  nH 


(2.47) 


for  the  components  of  the  moment  of  inertia  tensor,  /,  as  seen  in  the  reference  frame  of 
the  accelerometer  housing,  and  made  the  substitution 


^  V  ' 


(2.48) 


for  the  components  of  the  angular  velocity  of  the  housing,  as  seen  in  the  reference  frame 
of  the  rotating  housing. 

The  generalized  angular  forces  can  similarly  be  computed; 


N, 


I  p  _ 


dc 


F 

tj 


dRr, 


dR 


rtf  H 

p  9 


dV(R^,r^) 


•j 


dR 


+  \mRmj  +  A;, /z/j  + 


rl/[' 


kq 


OR 


Ou)'^ 


dR 


jt  (  ,H  .  ,P\ 

Ilq[^q  +  W  9  j  - 


1} 


dV{R^.r‘ 

dR^ 


+  2A„„  -f  , 


(2.49) 


where  we  have  used  the  symmetry  of  A,  and  made  the  same  substitution  for  the  conipo 
nents  of  I  and  u;^. 
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The  partial  derivatives  of  can  he  computed  hy  recognizing  (from  Ef}.  (2.31)) 


that  is  the  dual  of  ' ,  i.e., 


iP  _  *  qP  pP 
~  Ipm  "  tnn  ^‘pn- 


(2..30) 


The  partial  derivatives  can  now  be  computed 

1 


OR 


2  (Itm  Rm  j  ’ 


du 


,P 


OR 


bP 


.y  ^ttrn  R- 


nij- 


(2.51) 

(2.52) 


Substituting  these  expressions  into  the  equations  for  N  and  P  gives 


A',''’  .  5,... sr,,/;, (-■;'  +  *■'(•)- 


i)R 


■E  w'm  f riljSJ, , 
1 
2 


r:.’’  =  -  E<... 


q 


)■ 


(2.53) 


2  .^.1 


The  angular  equations  of  motion  ran  tiow  be  computed.  Setting  the  time  deriva 
tive  of  tlie  generaliz('d  mometita  equal  to  the  generalized  force  gi\('s 


.(lm,Prr,j  Uq  V  +  ^  1  j  - 


OV(R'\r’') 


,11  ,  ,P\ 


or: 


+  2A,„, 


+  2*  li’n^mji'lq  {^'q'  + '^'q  )  + 


Lagrange’s  equations  can  be  sitnplified  by  multiplying  both  sides  by  R'  and 
tuting  Eq.  (2.31)  into  the  result.  We  fitid 


(2.55) 

siibsti 


1  = 

u  ^Iq  \ 

'  ,H 

luJ 

+  2^‘<J^'lq  (^'q" 

+  u.') 

m  (R'\r'' 

or! 


4  2A,,  -f 


■mj 


This  ran  be  further  simplified  by  making  the  substitution; 

I'  —  At*,  t'  -i.  ,  I'  At' 

‘Iq  ~  '^in'lmn'nq  e  'qmn‘h,'^,n  • 


(2.5fi) 


(2.57 


and  applying  the  identity: 


We  find 


'f+o 


r  ,Pj,  (  ,H  .  ,P\  ,  iP  ,>  (  lH  ,  lP\  o 

I  hq  q  ^  q  )  i  ^jq  q  +  9  j  “ 


dV(R^,r’ 

orL 


_L  "  ,  I'  ,'  ‘  .'P  ^  I'  .  .'P  ,'P 

'  .^'^1  ‘qq'^q  ^tq^q'  ^Ic^cqm^m '*^q  ’ 


iP  l>  J 


fP  jH 


(2.58) 


+  2  A , , 


(2.59) 


This  set  of  nine  equations  can  be  viewed  as  a  .'{.xd  matrix  of  equations.  These 
equations  contain  the  six  forc<“s  of  constraint,  A,j,  that  force  the  components  of  to 
satisfy  the  condition  that  R^  be  an  orthogonal  transformation.  They  ran  be  eliminated 
by  separating  the  symmetric  and  antisymmetric  portions  of  our  matrix  of  equations. 
This  will  leave  us  with  six  equations  which  contain  the  constraints,  and  three  equations 
which  do  not.  Wo  find,  for  the  symmetric  part  of  our  matrix  of  equations. 


2A. 


V, . -n -.3  (-■'!' i;., (-'('+-’0 


+  - 
•) 


H 


dViR^’.r’’) 


;u 


—  +  R[ 


OR 


0\  {R'.r‘  ) 
0R\i 


(2.G0) 


These  equations  give  us  the  force's  of  constraint. 

Substituting  these  forces  back  into  our  nine  equations  of  iiiotion  eliminatc's  the 
symmetric  part  of  tlu’  matrix  of  equations,  raking  tlu'  dual  of  our  matrix  of  equations 
leaves  us  with  the  antisymim’trir  portion: 


fk, 


•>  /'  '  I' 

^  -'to-',  ^v-'v  -  ’k  '  ’ 
„  0V{R'’.r'') 

~  'klj  Rjy, 

OR  . 


,P  ,11 


(•2.G1 


This  can  be  simplified  by  usine  the  identity  in  F.q.  (2.58),  and  droi'ping  all  symmetric 
terms.  V\'e  find 


‘kq'^q  'kq'^  (, 


/II 

q 


/f  If  f  ,  wf 
'  .  /...«*'■  n  ^  J 


0\'[R'\r'’) 


or!’ 


'  2.G2 ) 


These  are  Euler’s  equations  in  the  reference  frame  of  the  accelerometer  housinp;. 

The  last  term  on  the  right  hand  side  of  this  equation  is  tlm  applied  toKjue  from 
the  coils  surrounding  the  proof  mass,  i.e., 


dV{R^,r^’} 

n  =  - 


dR 


To  see  this  more  clearly,  we  will  reduce  the  problem  to  two  dimensions,  and  forc''  all 
rotations  to  take  place  in  the  x-y  plane.  We  can  then  define 


R 


p  _ 


nj' 


(•2.r»i  I 


where  is  the  generator  of  rotations  about  the  r  axis,  and  has  roiii|)on<’rit.s  i. 

In  this  two-dimensional  sinqtlification,  \'{R^.r^  )  can  only  be  a  function  of  ft.  I  In' 
torque  generated  by  the  potential  T  is 


<)V 

7‘.;  = - .  Ci.ba) 

()ft 

We  can  rewrite  this,  using  tlie  rhain  rule  of  calculus: 


ov 

OV 

oRi:^ 

00 

ORi] 

Oft 

ov 

/  '  r’’  - 

hvi  mj 

/> 

-- 

orL 

ORi] 

Thus  Euler’s  equation  can  be  rewritten; 


j!  JP  _  !' 

-  ~‘kq 


‘1  _L  I'  -'P  -'P 

I  ‘  jq'^  q  '  ‘kq'q<;~^\  j 


n- 


■.  y’.tT 


The  components  of  the  moment  of  inertia  tensor.  I.  am  known  <'X])lir|' | v,  onl\ 
in  the  reference  frame  of  the  proof  m;iss.  Since  we  have  cho'-en  our  a\i>s  to  lie  alone,  t  hr 
principal  axes  of  the  proof  mass.  I  is  diagonal  in  this  refi'rem  ('  frann'.  If  t  In'  [>roof  mass  of 
the  accelerometer  is  comjilelely  symmetrical  so  that  all  throe  prim  ijral  moments  of  inert  i,i 
are  equal,  then  the  moment  of  inertia  tensor  will  be  diagonal  in  all  rotated  (’arte.-ian 
coordinate  systems.  I  bis  allows  us  to  simplify  the  right  hand  side  of  Euler's  eijuations 
If  the  proof  mass  was  not  completely  symmetrical,  then  term'  coupling  the  s(|uari'  ol 
the  angular  momentum,  of  the  housing  would  introdurc'  I'rror'  in  tlo'  measurement  of  the 
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angular  acceleration  of  the  accelerometer  housing.  Setting  1'^  =  IS,j  in  Eq.  (2.67),  we 
have 

luf  -  .  (2.68) 

This  equation  forms  the  cornerstone  for  an  understanding  of  the  detection  of  the  angular 
acceleration  of  the  housing  with  respect  to  inertial  space. 


2.2.3  Accelerometer  Equations  of  Motion 


We  have  now  completed  our  description  of  the  dynamics  of  the  SSA.  This 
description  does  not  yet  describe  how  a  measurement  of  the  orientation  of  the  proof 
mass  will  allow  us  to  measure  the  angular  acceleration  of  the  accelerometer  housing.  In 
order  to  complete  this  description,  it  is  necessary  to  pick  a  particular  set  of  angles  with 
which  to  describe  the  orientation  of  the  proof  mass.  This  will  then  allow  us  to  combine 
the  kinematics  of  this  description  with  Euler’s  equations,  and  derive  explicit  angular 
equations  of  motion  which  will  directly  relate  the  orientation  of  the  proof  mass  to  the 
angular  acceleration  of  the  accelerometer  housing. 

The  detection  of  linear  acceleration  is,  on  the  other  hand,  completely  described 
by  Eq.  (2.41): 


1  dV{R 

+  —  ' 
m 


P  .P\ 


=  al 


(2.69) 


where 


-  V$/r  -  x  x  X  x  r^.  (2.70) 

The  superconducting  coils  surrounding  the  proof  mass  generate  the  potential 
V.  The  position  of  the  proof  mass  with  respect  to  the  housing  is  controlled  by  these 
coils.  The  coils  are  rigidly  mounted  to  the  housing  of  the  accelerometer.  Ideally,  in  an 
accelerometer,  we  demand  that  this  applied,  i.e.,  measured,  acceleration  is  equal  ami 
opposite  to  the  acceleration  of  the  housing.  It  is  clear  from  Eq.  (2.70)  that,  in  order  to 
equate  the  linear  acceleration  of  the  housing  to  the  linear  acceleration  of  the  proof  ma.ss. 
we  must  keep  the  displacement  of  the  proof  mass  from  the  center  of  the  housing,  r^\ 
and  the  velocity  of  this  displacement,  r^,  to  a  minimum. 

This  can  be  achieved  by  using  force  rebalance  feedback.  This  type  of  feedback 
effectively  stiffens  the  magnetic  spring  between  the  housing  and  the  proof  mass,  and 
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at  the  same  time  dampens  any  oscillations.  The  applied  feedback  exactly  cancels  the 
acceleration  of  the  housing  and  holds  the  proof  mass  at  rest  with  respect  to  the  housing. 
This  forces  and  to  vanish. 

Although  Eq.  (2.68)  completely  describes  the  angular  motion  of  the  accelerome¬ 
ter,  it  is  not  clear  how  a  measurement  of  the  angular  position  of  the  proof  mass  describes 
the  angular  acceleration  of  the  housing.  The  orientation  of  the  proof  mass  of  the  SSA 
can  be  described  by  the  rotation  matrix  ,  but  such  a  description  is  not  convenient.  In 
order  to  explicitly  describe  the  motion  of  the  proof  mass,  it  is  necessary  to  parameterize 
the  rotation  matrix  in  terms  of  three  angles.  The  usual  choice  is  to  describe  R^  in  terms 
of  the  Euler  angles.  If  we  were  to  describe  R^  in  terms  of  Euler  angles,  then  we  could 
write  [5] 

R^  =  (2.71) 

where  J‘,  i  —  1,2,. 3,  are  the  generators  of  rotation  about  the  i-tli  axis,  and  have  com¬ 
ponents  J'^i^  =  —tjki-  The  Euler  parametrization  of  the  rotation  group  has  many  ad¬ 
vantages,  but  in  examining  Eq.  (2.71),  it  is  immediately  apparent  that,  for  small  angles, 
there  is  no  simple  relationship  between  rotation  about  the  x  axis,  and  0,  0,  and  0.  Thus, 
we  will  instead  choose  the  parameterization: 

rP  (2.72) 

This  places  rotations  about  all  three  axes  on  roughly  equal  footing.  Although  the  angles 
9x,  dy,  and  0^  do  not  exactly  correspond  to  rotations  about  the  x,  y,  and  z  axes  of  the 
accelerometer,  for  small  angles  this  is  a  good  approximation,  and  it  makes  the  succeeding 
calculations  more  tractable.  The  fact  that  the  three  successive  rotations  do  not  commute 
will  lead  to  equations  with  nonsymmetrical  cross  terms  in  the  second  and  higher  orders. 
However,  it  will  be  seen  that  these  unsymmetrical  terms  eventually  drop  out.  Expanding 
Eq.  (2.72)  into  explicit  matrix  form,  we  find 

cos0ycos02:  —cos0ysin9^  sinSy 

nP 

R  =  cos0isin0i  +  sin0isin0ycos0^  cos0^cos0.  -  sin0rsin0ysin0~  -sin0jrcos0y 
sin0isin02  —  cos0rsin0ycos02  cos0j-sin0yMn0,  +  sin0j-cos0,  cos0xcos0y 

(2.73) 
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iP 

U2  =  cos 


a  a  a 

Ox  — ~  —  sinWxCOS^u  -  , 

dt  ^dt 

tP  n  dO  ^  *  dOy 

Un  =  COSO'xCOSWu  — —  +SinPx  - -  ■ 

dt  dt 


(2.75) 

(2.76) 


These  equations  relate  the  rate  of  change  of  our  three  angular  coordinates  to  the  angular 
velocity, 

In  order  to  write  Euler’s  equations  explicitly  in  terms  of  ^x,  and  0^,  the 
torque  (Eq.  (2.63))  applied  by  the  superconducting  coils  surrounding  the  proof  mass 
must  be  related  to  partial  derivatives  of  the  potential  V  with  respect  to  0.  Oi.ce  more 
utilizing  the  chain  rule  of  calculus,  we  may  write 


dV 

dOx 


dV  dR^j 


dV 


d0x 


~  ^kml 


dv 


^  r> 


p 


(2.77) 


—  7’' 

—  ^  1  ’ 

where  we  have  substituted  our  parametrization  for  .  The  partial  derivative  of  the 
potential,  V ,  with  respect  to  0y  and  dz  can  similarly  be  computed: 
dv  dV  dR^j  dV 


dB. 


SzJ'  j2^-9zJ' 


dRlj  d9y 

dV  ,  ,  X 

(j^  cosBx  +  s\n0x)  R 

\  f  km 


R 


km 


m) 


dR], 

-  T2  cos  Bx  +  sin  ^x  1 


P 

mj 


(2.78) 


dv 

OBz 


dV  dR^^ 


dR^^  dB, 
dV 


dR^: 

dV 

dR^k, 

dv 


^9zj\eyj'^  ^-SzJ' 


JzJ'  (  rl  c; 


R 


km 


mj 


sin  By  +  cos  By^  e 


R 


km 


mj 


(j^  sin  By  -  sin  Bx  cos  By  +  J  ’  cos  Bx  cos  By)  R^ 

x  '  km  ■' 

sin  ByT'i  -  sin  0x  cos  ByT^  +  cos  Bx  cos  By'I'^.  ( 2.79) 


Solving  Eqs.  (2.77),  (2.78),  and  (2.79)  for  T[,  Tj,  and  T3,  we  find 

dV 


T[  = 


dBx 


(2.80) 
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dV 

II 

—  cos  9x 

d9y 

T'  — 

cos  9x 

dV 

^3  - 

cos  9y 

d9x 

c'e/j  “y 

.  ^  cos  6x  sin  9y 

sin 


CU5  Uy 

dV 

^  '  L.U&  Uj;  &lil  Uy  '■'  ' 

d9y  89, 

This  completes  the  parametrization  of  in  terms  of  9j.,  9y,  and 

We  can  now  compute  Euler’s  equations  in  terms  of  these  angles.  Substitutii 
the  expressions  for  T-  and  into  Eq.  (2.68),  and  solving  for  9j,,  9y,  and  9^,  we  find 


9x  tan  9yU)'^  -  sin  tan  9yu!^  — 

/  •  •  •  \ 

+  - ^  sin  9j;92  -  cos  9y  cos  9x9^  +  sin  9^  sin  9y9x  1 

cos  9y  ^  ' 

+  — ^  I  cos +  sin  9x  cos  9y9y  +  cos9x  sin  9y9x) - -^9y  sin  9y 

_  ^  ^  My  Sin  1  ^  ^  _  1  1  ^  ^  ^2.83) 

cos  9y  cos  9y  /  cos  9y  I  cos^  9y  ' 


9y  =  -  cos  9xUj''^  -  sin  9iU}'^  +  (u'^  cos  9^  -  sin  9x^  [9^  +  sin  9y9^ 


\  dV 

cos  9y92  -  -  -  +  9x9 2  cos  9y, 

I  d9y 


(2.84) 


e. 


cos  9x  .  ,ff  ,  sin0r  . 

(jJn  T  ^ 

COS  6y  cos  0y 


-  (uig" sin 9x  +  cos 9x)  — ^  f^r  +  sin 9y9^)  + 

^  >  cos  9..  \  I  cos  9„ 


cos  9y 

1  dV  dV 

7 - 20 - 

I  d9x 


cos  9, 


-9x9y. 


(2.85) 


Unlike  the  equations  for  the  linear  acceleration  of  the  six-axis  accelerometer, 
these  equations  which  relate  the  acceleration  of  9x,  9y,  and  9~  to  the  angular  acceleration 
of  the  accelerometer  housing  are  highly  nonlinear.  This  is  not  a  difficulty,  since  the 
proof  mass  of  the  accelerometer  is  surrounded  by  superconducting  coils  which  confine  its 
motion  and  force  all  9  and  9  terms  to  be  small.  Thus  we  can  linearize  our  equations  of 
motion,  and  drop  all  second  order  terms.  Expanding  Eqs.  (2.83),  (2.84),  and  (2.85)  to 
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first  order,  we  find 


1  dV 

^x+7  - 

^  de^ 
1  dV 
«v+7  — 

^  dffy 
1 

0^  +  j  — 
^  80, 


0  dV 

,•/«  4.  a  ,  ,  jHa  ,  ^  _ 

I  ^3  vy  “f*  u.^2  ^2  »  f  1 

0  dV 

-  +  -^ - . 


(2.86) 

(2.87) 

(2.88) 


These  equations  are  much  more  reasonable,  and  clearly  show  the  connection  between  the 
dynamics  of  the  proof  mass  and  the  angular  acceleration  of  the  housing.  It  is  clear  that, 
in  order  to  equate  the  angular  acceleration  of  the  housing  to  the  angular  acceleration  of 
the  angular  coordinates  of  the  proof  mass,  we  must  keep  the  angular  displacement  of  the 
proof  mass  from  the  center  of  the  housing,  0,  and  the  velocity  of  this  displacement,  0. 
to  a  minimum. 

This  can  be  achieved  by  applying  force  rebalance  feedback  to  the  angular  de¬ 
grees  of  freedom.  This  effectively  stiffens  the  dynamics  between  the  housing  and  the 
proof  mass,  and  at  the  same  time  dampens  any  oscillations.  The  applied  feedback  ex¬ 
actly  cancels  the  acceleration  of  the  housing  and  holds  the  proof  mass  at  rest  with  respect 
to  the  housing.  This  forces  6  and  $  to  vanish. 


2.3  Superconducting  Circuitry 

The  superconducting  inductance  bridge  has  been  analyzed  in  detail  by  Paik  [12] 
as  a  readout  circuit  for  a  gravitational- wave  transducer.  Input  and  output  impedances, 
and  forward  and  reverse  transductances  of  the  inductajice-bridge  transducer  have  been 
computed.  Here  we  carry  out  a  similar  analysis,  with  less  generality,  for  a  more  complex 
system:  a  six-axis  accelerometer  with  a  levitated  proof  mass. 

We  have  taken  our  analysis  as  far  as  we  are  able  without  restricting  the  design  to 
a  particular  geometric  configuration.  In  this  section  we  will  present  the  actual  geometry 
of  the  proof  mass  and  describe  the  positions  and  orientations  of  the  superconducting  coils 
which  surround  it.  We  will  also  examine  the  superconducting  circuitry  which  controls, 
levitates,  and  senses  the  position  and  orientation  of  the  proof  mass.  We  will  compute  the 
potential  energy,  V ,  of  the  superconducting  circuitry,  and  derive  the  transfer  functions 
that  relate  the  output  current  of  each  bridge  circuit  to  the  applied  displacement  of  the 
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proof  mass.  We  will  not,  however,  describe  any  construction  details,  but  rather  defer 
that  discussion  until  a  later  chapter. 

Figure  2.4  shows  a  perspective  view  of  the  proof  ma.ss  for  the  six-axis  accelerom¬ 
eter.  Each  of  the  three  intersecting  planes  which  make  up  the  proof  mass  is  defined  to 
have  a  thickness  of  2d,  and  a  length  of  2L.  Twenty  four  sensing  coils  and  twenty  four 
levitation  coils  are  mounted  on  eight  titanium  coil  forms,  as  shown  in  Fig.  2.4.  Each  of 
the  coil  form  surfaces  on  the  titanium  cubes  contain  one  sensing  coil  and  one  levitation 
coil.  Each  of  the  coils  in  the  accelerometer  is  labeled  according  to  its  position  with  re¬ 
spect  to  the  proof  mass  coordinate  system.  For  example,  is  the  inductance  of  the 
sensing  coil  above  the  xy  plane  in  the  -fi,  +y  quadrant  on  the  —z  side  of  the  proof  mass. 
Similarly,  is  the  inductance  of  the  levitation  coil  above  the  yz  plane  in  the  -y,  +z 
quadrant  on  the  -|-x  side  of  the  proof  mass. 

Each  of  the  coils  which  surrounds  the  proof  mass  remains  fixed  to  the  ac¬ 
celerometer  housing.  The  distance  from  the  center  of  each  coil  to  the  proof  mass  surface, 
and  the  orientation  of  each  coil  with  respect  to  the  proof  mass  surface  changes  as  the 
position  and  orientation  of  the  proof  mass  change.  This  forces  the  inductance  of  each 
coil  to  change  since  the  Meissner  effect  excludes  aU  magnetic  flux  from  the  body  of  the 
proof  mass.  Before  we  can  understand  how  the  superconducting  circuitry  functions,  it 
is  necessary  to  compute  the  inductance  of  each  superconducting  coil  in  the  SSA. 

Paik  [10]  has  shown  that  the  inductance  of  a  spiral  coil  above  a  superconducting 

plane  is 


L  =  Lo  +  Ar  =  A(ro  +  r). 


(2.89) 


where  A  =  A  is  the  area  of  the  spiral  coil,  tq  is  the  effective  spacing  of  the  coil 

above  the  superconducting  plane,  and  Ur  is  the  number  of  turns  per  unit  length  along 
the  radius.  This  expression  can  be  expanded  to  include  second  order  terms,  i.e., 

L  =  Lo  +  Ar-^r^  -^9\  (2.90) 

where  7  and  /3  are  nonlinear  coefficients  which  describe  the  change  in  inductance  through 
second  order  for  displacement  (r)  and  orientation  (0).  These  coefficents  will  play  a  major 
role  in  determining  the  resonance  frequencies  of  the  SSA.  These  coefficients  are  important 
when  the  size  of  the  levitation  or  sensing  coil  begins  to  approach  the  equilibrium  distance 
away  from  the  superconducting  plane.  We  will  calculate  those  parameters  in  Chapter  4 
for  our  particular  design. 


26 


Figure  2.4.  Perspective  view  of  the  six-skxis  accelerometer  proof  mass  and  coil  forms. 
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Table  2.1.  Initial  position  of  the  center  of  the  superconducting  coils  with 
respect  to  the  center  of  the  proof  mass  in  the  six-axis  accelerometer. 

Knowing  the  initial  position  and  orientation  of  each  coil  allows  us  to  compute 
the  new  spacing  and  orientation  of  each  coil  with  respect  to  the  proof  mass.  This  new 
spacing  and  orientation,  when  combined  with  Eq.  (2.90),  describes  the  inductance  of 
each  coil  in  the  SSA  as  a  function  of  the  position  and  orientation  of  the  proof  mass.  If 
x'  is  the  initial  position  of  the  center  of  a  coil,  then,  from  Eq.  (2.21),  the  coordinates  of 
that  point,  as  seen  in  the  cooi  dir.ale  system  attached  to  the  proof  mass,  arc 

X,  =  iZj-  (x'  -  rf)  .  (2.91) 

Only  one  of  these  three  coordinates  will  have  any  significance.  For  example,  in  order  to 
compute  the  inductance  of  the  L^yz  coil,  we  need  to  know  the  displacement  of  the  coil 
from  the  surface  of  the  proof  m<iss.  After  our  transformation  of  the  initial  position  of 
the  Ljzyz  coil,  only  the  z  component  will  be  needed,  as  the  new  position  in  x  and  y  will 
not  affect  the  inductance  of  the  coil,  since  the  Lxyz  coil  lies  parallel  to  the  x-y  plane. 

The  initial  position  of  the  center  of  each  of  the  coils  in  the  accelerometer  is 
given  in  Table  2.1,  where  d  is  the  initial  separation  of  each  coil  from  the  center  of  mass 
of  the  proof  mass  along  the  normal  direction  of  the  coil,  and  c  is  the  initial  separation 
of  each  coil  from  the  center  of  mass  along  the  direction  parallel  to  the  coil  surface. 

By  taking  the  inverse  of  Eq.  (2.73)  and  applying  Eq.  (2.91 )  to  the  initial  position 
of  the  center  of  each  coil,  and  expanding  through  second  order  in  9  and  r^,  it  is  possible 
to  calculate  the  displacement  of  each  coil  away  from  the  proof  mass  surface.  A  similar 
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calculation  yields  the  orientation  of  the  coil  with  respect  to  the  proof  mass.  Substituting 
the  displacement  and  orientation  into  Eq.  (2.90)  for  each  coil  in  the  SSA,  we  find 

Lxzy  —  Lq  +  A  "I"  c(^x  ~  ^z)  cByffz  ~  Tz^x  +  ^ x6 z  ""  "b  ^r)/^] 

-!['•«  +  -  6zf  -  2ryc{ex  -  0^)]  -  ^  (^x  +  ^l)  ,  (2.02) 

Lxzy  =  Lq  +  A^-Ty  +  ci-6x  -  ^z)  -  ce^6y  -  -  d{dl  +  0'l)/2^^ 

-  I  K  +  +  0x)]  -  ^  («x  +  ,  (2.9:3) 

Lxzy  =  Lq  +  A^-Vy  +  c{-0x  +  0z)  -  c6^dy  -  r::ej:  +  rx9z  -  d.(dl  +  9])/2]^ 

-  I  [^y  +  -  9,)^  +  2ryc(9,  -  0,)]  -^(9i  +  9^)  ,  (2.94) 

Lxzy  —  Lq  A  —Ty  4-  c{9j:  -|-  9^)  +  c9^9y  —  r^9j:  +  —  d{9^  +  ^?)/2 

-11^1  +  ^'(^x  +  Bzf  -  2ryc{9,  +  ^.)]  -  f  (^x  +  ^z)  <  (2.9.5) 

Lxzy  —  Lq  A  ^ry  4  c{—9x  4"  9^)  —  c9;9y  +  T^Ox  —  T^x^z  ~  f((^r  +  ^z)/2 

-  ^  H  -Ozf-  2ryci9x  -  ^z)]  -  ^  +  9l)  ,  (2.96j 

^zzy  =  +  A  ^fy  4-  c{9x  4-  6^z)  +  C^z^y  4  —  Tx^z  “  rf(^x  "b  ^z)/2| 

-  I  [^y  +  ‘^^(^x  +  ^z)^  +  2ryc(0x  +  ^^)]  -  f  (^X  +  ^z)  ,  (2.97) 

Lxzy  =  /-O  +  A  Ty  +  c(0x  -  ^z)  +  C^z^y  +  '■z(^i  -  ^x^z  -  <f(^r  +  ^z)/2] 

-  ^  [^y  +  c\9x  -  9xf  +  2ryC(9x  -  ^z)]  "  f  (^x  +  ^z)  ,  (2-98) 

Lxzy  =  Lo  +  A  Ty  +  c{-9x  -  9x)  ~  c9z9y  +  rx9x  -  rx9x  -  d{9l  +  9l)/2^ 

-  I  ['•y  +  c'(^x  +  0z)'  -  2ryc(^x  +  ~  I 

.^-xyz  =  f/o  +  A  -r,  +  c(^y  -  9x}  -  rx9y  +  Ty9x  -  d(6^x  +  ^y)/2] 

-  \  [r^z  +  c'(9x  -  9yf  4-  2r,c(^x  "  ^v)]  “  f  (^x  +  ^y)  •  (2.100) 

Lxyz  =  Lo  +  A  [-r,  +  r(-^y  -  ^x)  -  Tx^y  4  Ty^x  -  rf(^r  +  ^y)/2] 

-  I  ['•z  +  cH9x  4  ^y)''  4  2r,c(^x  +  ^y)]  "  f  (^x  4  ^y')  ,  (2.101) 

L^z  =  Lq  +  A^-Tx  +  c{-9y  +  9x)  -  rx9y  +  ry9x  -  d{9l  +  9l)/2 

-  I  ['•z  4  c\9x  -  9y  f  -  2r,c(0x  -  ^y)]  -  ^  (^x  4  9l)  .  (2.102) 

^lyz  =  f-o  4  A  —  Tj  +  c(9y  4  ^x)  —  ^x^y  4  ry^x  ~  d{9l  +  ^y)/2j 
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+  +  6yf  -  2rA9r  +  ^y)]  -  ^  {el  +  02)  ,  (2.103) 

Liyz  —  io  +  A  4-  c(~^y  +  ^r)  +  ^r^y  “■  ^y^J-  ~  d{9l  +  02)/2j 

-  J  -  ^y)'  +  2r,c(0,  -  ^v)J  -  f  (^x  +  ^y)  .  (2-104) 

Lxyj  =  Lq  -f  A  -f  c(6y  +  0jr)  +  rx0y  —  VyB^  —  d(0l  +  ^y)/2j 

-  I  [rl  +  c\er  +  Oyf  +  2r,c(0,  +  0,)]  -  ^  (ol  +  Ol)  ,  (2.105) 

Lx^  =  Zq  +  A  +  c(0y  —  0x)  +  ^xBy  —  Ty0x  —  (i{0l  +  02)/2| 

-  I  +  ^2(0,  -  0,)2  -  2r,c(0,  -  ^y)]  -  f  +  ^y)  ,  (2.106) 

■K'xp'  =  Aq  +  A  +  c{-0y  -  0J.)  +  rAy  -  ^y^x  —  d(0l  +  0y)/2j 

-  I  (/z  +  c\0x  +  Oyf  -  2r.r(0,  +  0,)]  -  ^  (^x  +  ^y)  .  (2-107) 

iyzx  =  1^0  +  A  [-r,  +  c(0z  -  0^)  +  c(0s,0,  +  0,0J  -  r,0,  -f  r,0,  -  ^^02  +  02)/2] 

-11^1  +  +  2r,c(0,  -  ^z)]  -  ^  +  ^z')  .  (2-108) 

Lyzx  =  Lo  +  A  -Tx  -  c{0z  4-  0y)  4-  c{0z0x  -  0y0x)  -  ry0 z  4  rz0y  -  c((0y  +  0^/2] 

-11^1  +  +  2''xC(0y  +  0z)]  -  f  (^y  +  ^z)  .  (2-109) 

=  Io4-A[-r^-fc(0y-0z)'C(0y0x  +  0z0x)-rj,0z4r..0, -rf(02  +  02)/2 

-  I  [i-x  +  c'(^y  -  ^z)"  -  2rxC(0,  -  ^z)]  -  ^  {el  4  0^)  ,  (2.110) 

Lyzx  —  ^0  4-  A  [“I'x  4  c{0z  4  0y)  +  c{0y0x  —  0z0z)  —  I'y^z  +  I'z^y  ~  l^(^y  ^  ^z)/2j 

-  J  [i-x  +  i''(^y  -  ^z)"  -  27-,c(0,  +  (»z))  -  ^  4-  el)  ,  (2.111) 

^yzx  =  -^-0  4-  A  r.r  4  c(0y  -  0z)  -  c(0y0x  +  ^z^x)  +  l-y^z  -  I'z^y  -  <^(^y  +  ^z)/2] 

-  I  {rl  +  cHey  -  0,f  +  2rxC(0y  -  0^)]  -  ^  (»y  +  ^?)  ,  (2-112) 

Lyzx  =  Lq  ■\-  A  r X  -\-  c(0z  +  0y)  4  c{0y0i  —  Ozei)  +  I'y^z  ~  I'z^y  “  l^(^y  4  el)/2 

-  J  [i-x  +  i’'(^y  4  0z)'  4  2rxC(0y  4  0z)]  -  ^  {el  4  0?),  (2.113) 

/-y?!  =  io  4  A  T-x  4  c(0z  -  0y)  4  4^y^x  4  0z^x)  4  ^y0z  -  rz05^  -  (f(02  4-  02)/2 

-  J  [i-x  +  I^'(^y  -  ez?  -  2rxC(0y  -  ^z)]  -  ^  {el  4  02)  ,  (2.1 14) 

Lyjf  -  Lo  +  A  rx  -  r(0j  4- 0y)  +  f(0,0x  -  0y0x)  4  ry0z  -  rz0y  -  d(02  4  0?)/2] 


.'iO 


+  ^  (^y  +  ■■?)  ’  (-115) 

where  we  have  simplified  our  notation  from  to  r,  since,  from  this  point  on,  we  will  only 
be  interested  in  displacements  of  the  proof  mass  relative  to  the  accelerometer  housing. 
The  fact  that  the  equations  for  inductance  are  not  symmetric  stems  from  the  non 
commutivity  of  the  three  successive  rotations  used  to  describe  the  orientation  of  the 
proof  mass  in  Eq.  (2.72). 

Although  each  coil  form  in  the  accelerometer  contains  two  coils,  they  share 
a  common  center,  and  we  have  been  able  to  treat  them  identically.  Thus,  both  the 
levitation  and  the  sensing  coils  obey  the  same  formula  for  inductance  but  the  parameters 
Z/Q,  A,  7,  and  d  are  different,  so,  later  in  our  formulas,  we  will  distinguish  these  coils 
by  substituting  Ls,  As,  /3s,  7s,  and  ds  for  the  sensing  inductance  parameters,  and  L/,. 

I^Li  lLi  and  di  for  the  levitation  inductance  parameters. 

The  SSA  is  oriented  in  the  so  called  “umbrella”  orientation,  so  that  all  three 
sensitive  axes  make  the  same  angle  to  the  vertical.  This  orientation  is  identical  to  the 
orientation  of  the  gravity  gradiometer,  and  has  the  advantage  of  distributing  the  pull  of 
gravity  equally  among  the  three  axes. 

2.3.1  Levitation  Circuitry 

In  the  presence  of  gravity,  the  r^,  Ty,  and  r,  levitation  circuits  each  must  equally 
levitate  the  proof  mass.  These  three  circuits  arc  shown  in  f'ig.  2.5.  In  e.\amining  the  r^ 
levitation  circuit,  we  see  that  all  four  of  the  l‘'vitation  coils  are  located  at  the  back  of  the 
X  plane  of  the  proof  mass.  If  a  current  is  stored  in  the  levitation  circuit,  all  four  coils 
will  push  the  proof  mass  upward  equally.  This  is  the  source  of  the  force  that  levitates 
the  proof  mass.  The  additional  transformer  is  included  simply  so  that  we  can  add  or 
subtract  current  from  the  circuit  in  order  to  raise  or  lower  the  proof  mass  position  and 
apply  feedback. 

The  7-j,  levitation  circuit  contains  a  single  loop  about  which  the  trapped  flux 
must  remain  constant.  We  can  treat  the  five  inductances  in  this  circuit  as  one  equivalent 
inductance: 

f-yzr  +  f'yz3r  +  "h  ^"y^  +  (2.116) 

If  we  store  a  levitation  current,  in  this  circuit,  then  tlu'  amount  of  trapped  flux  is 
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(4Li  +  L)Ir,,  and  the  potential  energy  of  this  circuit  is,  therefore, 


(2.117) 


Substituting  the  required  inductances,  and  expanding  this  to  second  order  in  r  and  0 
gives 


V'r.  = 


+  L)ll 


‘ilrMTr  + 


8Ai 


ALi,  +  L 


+  7L 


+  ‘lI^ALiSyr,  -  0,ry)  +  +  c^lL  +  Hl)  •  (2.118) 


Thus  we  see  that  our  r^,  levitation  circuit  provides  a  levitation  force,  and  an 

additional  spring  constant,  2(8A^/(4Z./,  +  L)  +  7^,]/^^  ^'here  is  also  an  interesting  cross 
coupling  term. 

Similar  analysis  of  the  Vy  and  levitation  circuits  leads  to 


Vr..  = 


(‘^Ll  +  1)1} 


-  + 


SAl 


+  ir\r; 


2 

+  2i\AL(e,r,  -  e,T,)  +  (A/,rf,  +  r'u  +  i3i)  [el  +  e:) .  (2.1 19) 


Vr. 


8A 


2  ^  +  uj  l;j] 

+  2IIAl  {9,ry  -  dyv,)  +  [AlcIl  +  +  ^h)  (oj  +  O'^)  .  (2.120) 


The  three  angular  levitation  circuits  are  a  bit  more  complicated.  These  circuits 
are  shown  in  Fig.  2.6.  Each  of  the  three  circuits  contains  four  levitation  coils,  and  one 
feedback  transformer.  Each  set  of  four  levitation  coils  is  mounted  on  the  upper  side  of 
one  of  the  proof  mass  planes.  The  circuit  contains  the  four  levitation  coils  on  the 
upper  side  of  the  proof  mass’s  c  plane.  Two  independent  currents  arc  stored  in  each  of 
the  two  superconducting  loops  in  the  angular  circuits.  In  the  circuit  shown  in  Fig. 
2. .5,  the  current  in  the  loop  on  the  left,  tends  to  rotate  the  proof  lua.ss  about  r  axis 
in  a  counter-clockwise  manner,  and  the  current  in  the  loop  on  the  right.  i2e.~  tends  to 
rotate  the  proof  mass  about  the  x  axis  in  a  clockwise  manner.  Both  currents  tend  to 
push  the  proof  mass  downward  in  the  — c  direction.  I'he  feedback  transformer  is  included 
so  that  the  angular  position  of  the  proof  mass  can  be  controlled.  If  the  same  current  is 
stored  in  each  loop  of  the  circuit,  then  any  additional  current  added  by  the  feedback 
transformer  will  tend  to  decrease  one  current  and  increase  the  c-iher.  I  his  change  in 
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the  currents  in  the  two  loops  will  apply  a  torque  to  the  proof  mass.  This  ran  Ik  seer, 
quantitatively  by  computing  the  potential  energy  of  the  Oj.  levitation  circuit. 

The  two  flux  conser.ation  equations  for  the  6j.  circuit  are.  ignoring 

ihOz  +  h)(Lzyz  +  J-'ryz)  +  +  *1  —  =  2/./,/)^^  +  (ho^  —  12.121  I 


{hei  +  h)i^xyz  +  ^^ryz)  +  +  *2  “  ^Ifii  “  —  2A/,/2^j  +  ihSj  “  (2.122  i 

The  potential  energy  of  the  levitation  circuit  is  simply  the  olectromagnetic  eneig> 
stored  in  the  inductances: 

''fii  “  -j[l-ryz  +  )(-fl(?x  +  T  +  l-ryzKho^  +  '2)‘' 

+  2^Ai26^  +  ^2  ~  “  T)^-  (2.12.1/ 

Solving  Eqs.  (2.121 )  and  (2.122)  for  r'l  and  i2,  substituting  these  current.s  and  the  required 
inductances  into  our  expression  for  and  expanding  to  second  order  in  r  and  6 
gives 


+  hej)^L  +  2^126^  -  hoi)  I- 
+  (^20 r  +  +  (fl?,  “ 


+  (^29x  +  A^Sx  ■ 
+  ( ^26^  ~  ^le,:  . 


Ll  +  L  ^  ‘ 

'\i(2LL  +  L) 

1‘iALi,  +  1) 


'■j  +  ( A«i  +  A^x 


2Li(Li  +  L) 


+  ")/-  +  (  A^^’r  ^  2i>  A  i[(i  yr  j.  —  9^ry] 


+  {4x  +  +  (A.tx  - 

+  Ufk  +  A^^x )  (a/tA.  +  c'li.  +  %)  {o'i  +  9l)  ■ 


(2.121) 


This  potential  looks  (piite  complicated;  however,  wi-  can  iinim dia'ely  recognize  the  main 
features.  The  difference  in  the  two  -iirrents  /jo^  and  129^  ap()lies  a  torque  to  the  prooi 
mass  about  the  x  axis,  and  both  currents  push  the  proof  mass  downward  in  the  -r 
direction.  The  rest  of  the  terms  describe  the  “magnetic’'  spring  du<'  to  the  circuitry. 

Similar  analysis  of  the  9y  atid  9;  levitation  circuits  gives 


+  (A^^iy  +  l{g^)Air^  +  -  /]A^)c.\/,^y 


+ 


^  Ll-\-L 


+  (hsy  + 


{^2By  -  J\By)  IJ  ^  {^\By  ^2By)^lX^z^y  -  ^y'^z) 

+  ('J..  +  +  ('2«,  -  ''>^'i‘2iaiX+Lfy 

+  2^^2By  +  ^iBy)  {^idi  +  C^lL  +  0l) 


il2By  -  lfBy)ALC0^e, 


Vb,  —  {JIb,  ^IB.)^L  +  -^{hB.  -  hBrf  ly 

+  +  ^lBz)^L'’'y  +  {I26,  -  ^xbM^L^z 

+  (ik  +  +  ('m.  +  i/t 


(2.125) 


+ (/2.,  +  'f.,  )7^«; + (/2..  -  /,». 

+  [^idt  +  c^7l  +  /3£,^  . 


(2.126) 


We  can  now  combine  all  the  different  potentials  applied  to  the  proof  mass  by 
the  superconducting  levitation  circuitry.  Because  the  SSA  is  operated  in  an  umbrella 
orientation,  equal  amounts  of  levitation  current  are  stored  in  the  three  linear  circuits, 
thus  we  will  let  Ir^  =  hy  -  Ir.  =  f+  in  Eqs.  (2.118),  (2.119),  and  (2.120).  Similarly, 
we  will  store  equal  amounts  of  levitation  current  in  each  of  the  three  angular  circuits. 
Since  we  do  not  initially  want  to  apply  any  torque  to  the  proof  mass,  we  will  store  equal 
amounts  of  current  on  each  side  of  the  angular  circuits.  With  this  in  mind,  we  will  let 
hB,  =  hB,  =  IxBy  =  hBy  =  hs,  =  hB,  =  /-  in  Eqs.  (2.124),  (2.125),  and  (2.126). 
Combining  all  these  equations  together,  we  find  that  the  total  potential  energy  of  all  the 
levitation  circuits  together,  Vi,  is 

Vl  =  + L)II+(HILl  +  2{II- ll)At[r,  +  ry  +  r,) 


+  Uni  +  +  c^lL  +  0l)  +  2ll  [01  +  91  +  9] 

+  2(ll  -  ll)\i^  -  Ty)  +  9y(r^  -  r,)  +  9 Ary  -  r^)] . 


(2.127) 


2.3.2  Sensing  Circuitry 


We  have  claimed  that  the  sensing  circuitry  in  Fig.  2.7  will  provide  us  with 
information  about  each  of  the  six  degrees  of  freedom.  At  this  point  we  will  analyze  the 
sensitivity  to  motion  in  the  degree  of  freedom.  The  sensing  circuit  can  be  simplified 
to  the  circuit  shown  in  Fig.  2.8,  where  we  have  replaced  the  other  sensing  inductance 
bridges  by  their  equivalent  inductance  at  equilibrium,  Ls-  The  load  on  the  portion  of 
the  circuit  sensitive  to  is  then  Lsq  +  5Ls,  where  Lsq  is  the  input  inductance  of  tlie 
SQUID  amplifier. 

The  sensing  circuit  for  r^.  in  this  simplified  form,  contains  two  independent 
superconducting  loops  about  which  the  trapped  flux  must  remain  constant.  In  tliis  Cdsc. 
the  initial  flux  is  zero,  and  we  may  write  two  flux  conservation  equations: 


{Lyzx  T  iyzxlf  T  ^^yzx  IoCOS{u)rzf)  +  (Lsq  +  =  0, 


(2.128) 


(Ty2x“t“f'y2x)f  Lyzx  "i"  f^y2x)^rj  T  L'yzx^Q  ros(u.v j f)  —  0*  (2.129) 

Solving  these  two  equations  for  the  current  through  the  SQUID  amplifier,  we  find 
^  _  _ /q  COs(u>fjt)  (Z>y2x.^V2r  ~  ■f'yzr  f'yzT  ) _ 

[(^SQ  "I"  ^Ls){Lyzx  T  I'yzx  T  ^jTjx  T  f'yTx)  "I”  {Lyzx  T  ^yixlff'yei  "f"  f'ylx)] 

(2.130) 

Substituting  Eqs.  (2.113),  (2.108),  (2.110),  and  (2.115)  for  Lyzi,  Lyzx<,  L^x^  and  Lyjx^ 
respectively,  into  the  above  equation,  we  find,  after  expanding  to  second  order  in  r  and 

0, 


irz  =  /ocos(u;x^t) 


As^x  Asl^ygx  -  TzOy)  _  C^(A|  +  JsLsWy^Z 

Lsq+6Ls  Tsq  +  6Z/5  Ls(Lsq  + 


(2.131; 


Thus  we  see  that  jVi  is  directly  proportional  to  the  position  of  the  accelerometer  in 
Similar  analysis  of  the  y,  z,  6j^,  9y,  and  sensing  circuits  leads  to 


i  =  Iocos{ur^t) 


iV,  =  /ocos(a;r,0 


ifl,  =  Iocos{ugj) 


Asl^z^x  -  r^dz)  r^(.\^’  +  lsl^s)Sx0z 


Lsq  +  6Ls  Lsq  +  6/is 


Ls(  Lsq  +  6A5) 


+  •■•  , 


(2.132) 


Asr-z  ,  As(rx9y  -  Tyg^)  r^(A|  +  isLsWxOy 


_ ^  _  I  »  y 

Lsq  +  6L5  Lsq  +  6L5 


Ls(Lsq  +  ^Ls) 


+  ...  , 


cAs^x  (A|  +  jsLsyOyrz 
Lsq  +  ^Ls  Ls(Lsq  +  GLs) 


(2.133) 

(2.1.34) 


Figure  2.7.  Supcrcon 
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Figure  2.8.  Simplified  superconducting  sensing  circuitry  for  the  degree 

(A^  "F  'y s X  x^z 


of  freedom. 

isy  =  Iocos{ugJ) 

cAsOy 

Lsq  +  6Z.5 

19,  =  ZoCOS(W9,<) 

cAsOz 

Lsq  +  6Ts 

Lsi^sQ  +  CZ.s) 
.  (A|  +  lsLs)c9xry 

•  r  /  T  .r>r\  r*** 


+  ... 


(2.13.5) 

(2.136) 


The  final  current  that  flows  through  the  SQUID  is  the  combined  total,  iV, 
+  tVy  +  if,  +  ie^  +  igj,  +  ie,,  and  if  each  bridge  is  modulated  with  a  different  frequency 
current,  then  it  is  possible  to  recover  each  individual  output.  This  is  accomplished  by 
utilizing  six  lock-in  amplifiers,  each  one  tuned  to  the  separate  carrier  frequency. 

Each  sensing  circuit  will  affect  the  dynamics  of  the  proof  mass.  As  the  proof 
mass  changes  position  with  respect  to  the  bridge  inductances,  the  currents  flowing 
through  the  circuit  must  redistribute.  This  change  in  current  patterns  redistributes 
the  energy  stored  in  the  sensing  bridge  and  changes  the  magnetic  potential  energy.  .\s 
the  amount  of  sensing  current  is  increased,  the  stiffness  of  the  “magnetic”  springs  sur¬ 
rounding  the  proof  mass  is  increased.  In  order  to  understand  this  quantitatively,  it  is 
necessary  to  compute  the  amount  of  energy  stored  in  the  sensing  bridges  as  the  position 
of  tlie  proof  mass  changes.  The  total  electromagnetic  energy  stored  in  the  sensing 
circuit  is,  from  Fig.  2.8, 

^  ~  2  2  ( Z  T  Zo  cos  u)tj  t)  F  ~  Ly^x  (Z  —  iV, ) 
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+  +  loCOSUJr^t  —  ir^)^  +  “  (l^SQ  +  5Z/s)  (2.137) 

This  is  not  the  correct  potential  energy  of  the  sensing  circuit,  however.  As  the 
proof  mass  moves,  the  current  source  attached  to  the  sensing  circuit  must  maintain  a 
current  /ocosu^t.  This  forces  the  current  source  to  supply  energy.  The  voltage  induced 
across  the  current  source  is  the  time  rate  of  change  of  the  magnetic  flux  across  the  bridge: 

^  ^  [Lyzxil  +  loCOSUrj)  +  Lyji{I  +  loCOSUrJ  -  )]  •  (2.138) 


The  power  supplied  by  the  current  source  at  any  instance  of  time  is  the  product  of  the 
voltage  drop  and  the  supplied  current,  i.e.. 


P  -fo  cos  [//y2j,(/ -j“  To  cos t )  "f"  ^yzxi^  "b  /q  COS  UJ,.^  t  ^rx  )]  * 


(2.139) 


If  we  now  displace  the  proof  mass  a  small  distance  6x  away  from  equilibrium,  then 
we  can  compute  the  force  in  the  x  direction,  fx-  The  energy  we  have  received  in  this 
displacement  is  fx6x.  The  change  in  the  total  electromagnetic  energy  stored  in  the 
sensing  circuit  is  {dEldx]6x,  and  the  amount  of  work  done  by  the  current  supply  is 


Wc 


=  I 

J"  -fo  cos  W^rx  ^  [  ^yzxiy  ^  T  ^0  COS  )  -j-  Z'yJx  (  ^  "b  ^0  COS  tt^rx  ^  Vx  )] 

^ /q  COS  Ui^Xx  ^  ^  [Eyzx{^  /q  COS  U^Xx  0  d"  Tyj^(/  "b  /o  COS  U^Xx  ^ 
d 

=  loZOSUrJ—  [Lyzx{I  -k-  loCOSi^rj)  +  Lyzx{I  +  loCOSUJrJ  “  Irx)]<5x.  (2.140) 


Thus  conservation  of  total  energy  demands  that 


Wc  =  -b  frSx. 

dx 


(2.141) 


Substituting  in  Eq.  (2.140)  for  Wc  and  solving  for  Z^,  we  find 


/x  =  -  loCOSU^rJ  [Lyzx{I  +  loCOSUrJ)  +  LyJxil  +  loCOSUrJ  -  Zxx)])  • 


(2.142) 


Thus  we  .see  that  the  proper  potential  energy  of  the  sensing  circuit  is 


V^rjt  ^  A)  COS  /  -f"  A)  cos  U^rx  0"^  ^yzxi,^  ^0  COS  ^  )]  • 


(2.143) 
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Combining  like  terms  and  simplifying,  we  find 

“■  2  ^yzx  L-^x  d"  ~  ^  i^yzx  4"  Lyjx)  Iq  cos^u^n^ 

4"  2  4"  t'yzx  4-  Lsq  4-  5X5)  4-  Ly^)  /iVx-  (2.144) 

Solving  Eqs.  ^2.128)  and  (2.12^)  for  f,  substituting  /,  and  Lqs. 

(2.108),  (2.110),  and  (2.115)  for  Xy^i,  X^i,  and  Xy^j,  respectively,  into  the  above 
equation,  we  find,  after  expanding  to  second  order  in  r  and  6, 

1 


Vr^  =  -/o(l  4- cos2a;r,«) 


-Ls  + 


( 


A| 


75 


Lsq  4-  6X5  ^  2  j  ^ 


jfj  +  \  i^sds  4-  cS5  4-  0s)  {el  4-  02) 


(2.145) 


This  shows  the  relationship  between  the  square  of  the  sensing  current  and  the  “stiffness” 
of  the  accelerometer.  The  sensing  current  oscillates  at  a  high  frequency,  well  above 
the  resonance  frequencies  of  the  accelerometer,  and  it  is  an  excellent  approximation  in 
practice  to  ignore  the  time  varying  portion  of  the  potential. 

Similar  analysis  of  the  fy,  r^,  0^.,  0y,  and  0^  gives 

1 


^ry  =  ^/o(l  4- cos2u;ry0 


-Ls  4- 


A5 


4- 


7s 


4- 


^A|' 


Ls 


Lsq  4-  6Xs  2  j  ^ 
e]  4-  ^  (A5d5  +  c'^is  +  0s)  {el  +  el) 


(2.146) 


Vr.  =  ^^0(14-  COS  2Wrx  <  ) 


-Ls  4- 


A5 


Lsq  4-  6X5  2 


4-?^  r2 


(2.147) 


=  ;^/o(l  4- cos2u;9^t) 


-Ls  4- 


c^Al 


Lsq  4-  6X5 


+  j  5  (Asrfs  +  c27s  +  0s)  {el  +  02)  + 


(2.148) 


=  ^^o(l  4- cos2w9^f) 


-Ls  4- 


4- 


'c2A|' 


Lsq  4-  6X5  j 
4-  -  (A5d5  4-  C275  +  /Is)  -I-  02^  q. 


(2.149) 


'^6.  =  -/o(l  4- cos2w«,<) 


-Xs4- 


Lsq  4-  6X5 
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+  (^)  \  +  T"v  •  (2.150) 

In  our  analysis  up  to  this  point,  we  have  driven  all  sensing  circuits  with  the  same 
magnitude  of  sensing  current,  Iq.  It  is  desirable  to  separate  the  effects  of  the  linear  and 
angular  sensing  circuits  by  driving  each  with  a  differrct  mcgcituHe  of  oscillating  "'urrent. 
Thus,  we  will  drive  all  three  linear  circuits  with  an  oscillating  current  of  magnitude  7^, 
and  all  three  angular  circuits  with  an  oscillating  current  of  magnitude  Ig.  The  total 
potential  derived  from  all  six  sensing  circuits.  Vs,  is  then 


+  +  irl  +  rl  +  rl) 


Vs  =  - -{1,^  +  Ig^)Ls  +  -  7/ 


4  1  Lsq  +  6Ls 


(2.151) 


where  we  have  dropped  the  time  dependent  portion  of  the  potential. 

The  complete  potential  for  the  SSA  is,  from  Eqs.  (2.127)  and  (2.151), 

V  =  VL  +  Vs  =  Vo  +  2{ll-ll)AL{rr  +  ry-^r,) 

+  2(11  -  /+)Al  -  Tj,)  +  Oyir^  -  rj  +  0,{ry  -  r^,)] 


+  2  llj^  +  ll 

+  1  [/2 - ^ 

“  >4  T  ,  a 


2 


(rl  +  rl  +  rl) 


4  [  Lsq  +  6Ls 


+  (4  +  )‘7r  {^x  V  r  +  r^) 


+  2(ll+ll)  [A^dL  +  C^IL  +  Hl)  +  2ll 


Ll  +  L 


+  4  ilr+Ie) 


+  Asds  +  c^7s  +  (is  + 


Lsq  +  6Ls 


(01  +  01  +  02), 

(2.152) 


where 


To  =  ^(47,£  +  L)ll+6llLi^  -  ^(72  +  7,2)7s. 


We  can  rewrite  the  potential  V,  by  defining 


(2.153) 


foe  =  2{II-II)Al, 


(2.154) 

(2.1.55) 
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ki  =  4 


Ti  =  4 


^-Ll  ^  +  [4Ll-\-L'^  2  jj 


(4  +  /!)  (Airfz,  +  c^U  +  Pl)  +  II 


r.  -  ^  {(/?  +  /!} 


Ll  +  L 

,  A  I  .2._  ,  /3,1 


(2.156) 


(2.157) 


A5_ _ \  ro2.-«> 

.  X  -  V.  y 


Ls  ‘  ■  •''  •  ""J  '  •'’isQ  +  eisJ 

We  find,  after  substituting  these  expressions  into  our  potential  V , 

V"  =  Vb  -  fDc{rx  +  Tj,  +  r^) 

-  foe  [^xi^z  -  Ty)  +  Oy(rx  -  r^)  +  0^{ry  -  r^.)] 

+  2^^^  ^5)(^X  +  +  2^^^  d-  Ts){^x  +  +  &l)-  (2.159) 

It  is  immediately  clear  how  each  portion  of  the  potential  contributes  to  ♦he 
dynamics  of  the  SSA:  k^,  and  kg  are  the  linear  spring  constants  generated  by  the  lev¬ 
itation  and  sensing  circuitry,  and  Tg  are  the  angular  spring  constants  generated  by 
the  levitation  and  sensing  circuitry,  foe  is  the  DC  force  generated  by  the  levitation 
circuitry,  and  foe  is  also  a  measure  of  the  cross  coupling  between  the  linear  and  angular 
levitation  circuits. 


2.4  Six- Axis  Accelerometer  Transfer  Functions 

In  the  last  two  sections,  we  have  derived  the  equations  of  motion  for  the  SSA, 
the  sensitivity  of  the  sensing  circuitry  to  changes  in  the  position  of  the  proof  mass,  and  the 
potential  energy  of  the  magnetic  springs  which  levitate  the  proof  mass.  In  this  section,  we 
will  combine  all  the  information  in  the  previous  sections  and  describe  the  SSA  in  terms 
of  transfer  functions  relating  acceleration  to  displacement  and  displacement  to  output 
current.  This  will  enable  us  to  describe  the  fundamental  noise  of  the  accelerometer  and 
the  potential  sensitivity  of  the  SSA  in  the  next  section. 

A  displacement  of  the  proof  mass  in  the  SSA  away  from  equibbrium  is  echoed 
by  a  current  response  i,  in  the  input  coil  of  the  SQUID  amplifier.  The  six  t;  nsfer 
functions  relating  the  output  current  of  the  SSA  to  the  position  and  orientation  of  the 
proof  mass  can  be  derived  from  Eqs.  (2.131)  through  (2.136).  We  find: 

Hr,,  =  —  =  IrCOSUrJ- - ,  (2.160) 

Tj.  Lgq  +  bLg 


43 


=  IrCOSUr^t—^^f-—,  (2.161) 

Ty  Lsq  +  6Ls 

//r,.  =  —  =  IrCOSUrJ- - (2.162) 

rz  Lsq  +  bLs 

Ilg^,  =  ^  Igcosuej  ^  ,  (2.163) 

Ox  Lsq  +  01/5 

^  ^  fecos^eyt---  ,  (2.164) 

//5Q  +  oZ/5 

^  ^  i^cosu>eJ—~/—,  (2.165) 

Oz  Lsq  +  61/5 

where  we  have  dropped  all  second  order  terms,  and  replaced  Iq  by  fr  and  Ig  in  the  linear 
and  angular  sensing  circuits,  respectively. 

An  acceleration  of  the  SSA  is  echoed  by  a  displacement  of  the  proof  mass.  This 
displacement  of  the  proof  mass  is  described  by  the  equations  of  motion  for  the  SSA, 
Eq.  (2.69)  for  the  linear  degrees  of  freedom,  and  Eqs.  (2.86),  (2.87),  and  (2.88)  for  the 
angular  degrees  of  freedom.  Substituting  the  potential,  V',  from  Eq.  (2.159)  into  these 
equations,  we  find: 


Ti+UJ^Tx  = 

^  -  «,)  +  «f . 

(2.166) 

Ty  +  U^ry  = 

^  +  ^(9x-9x)  +  a^, 
mm  “ 

(2.167) 

r'j  +  = 

—  +  — (0x-^v)  +  af, 
mm 

(2.168) 

9x  +  = 

j  (rz  Ty)  +  u;'/' +  9yu;'j^' 

-  wf 

(2.169) 

9y  +  (Jjg9y  = 

^(rx-rz)-l-d’f -9xd>','f 

—  u)'/^  9  z  +  9x, 

(2.170) 

9z  +  ‘^g9z  = 

^  (fy  rx)  +  +  9xU>2^ 

—  9x  +  9y, 

(2.171) 

where  we  have  dropped  the  superscript  from  r^,  dropped  all  second  order  terms,  and 
defined 


^1/  + 
m 

tl  4-  Ts 
I 


(2.172) 

(2.173) 


The  SSA  is  orientated  in  the  “umbrella"  orientation  so  that  all  three  axes  make 
the  same  angle  with  the  vertical,  and  the  acceleration  of  gravity  is  applied  equally  to 
the  three  axes.  With  the  SSA  in  the  “umbrella”  orientation  on  the  surface  of  the  earth. 
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levitation  currents  are  stored  in  the  SSA  so  that 


fpc  _  ^ 

m  v/3 

This  levitation  force  cancels  the  DC  component  of  the  external  acceleration, 
tions  (2.166'!  through  (2.171)  can  be  rewritten; 


(2.174) 
.  Eqna- 


+  =  a,,  (2. 17.^1 

+tL>g0,  =  o,,  (2.176) 


where 


flx 

= 

„e,9e  9e 

+  75  +  75* » ■ 

Oz), 

(2.177) 

dy 

= 

75  ^  75*^-  “ 

9x), 

(2.178) 

dz 

= 

„e,9e  9e  . 

+  75  75**'  ■ 

(2.179) 

Otx 

= 

+  OyU)'^^  —  U)^  l>y  + 

(2.180) 

= 

r,)fc.2 

- 

(2.181) 

= 

+  +  ^'"Oy. 

(2.182) 

Up  until  this  point,  we  have  assumed  that  no  damping  l'I  the  motion  of  the 
proof  mass  occurs.  In  reality,  the  superconducting  circuitry  and  the  e..^  irimental  support 
structure  force  the  Q  of  the  modes  in  the  SSA  to  be  finite.  The  damping  of  the  proof  mass 
in  the  linear  degrees  of  freedom  can  be  modeled  by  the  addition  of  a  velocity  dependent 
force,  ruirlQr  lo  the  equations  of  motion  for  the  SSA,  where  is  the  quality  of  the 
linear  resonance.  The  damping  of  the  proof  mass  in  the  angular  degrees  of  freedom  can  be 
similarly  modeled  by  the  addition  of  a  velocity  dependent  force,  dijgIQg  to  the  angular 
equations  of  motion  for  the  SSA,  where  Qe  is  the  quality  of  the  angular  resonances. 
Adding  these  additional  terms  to  Eqs.  (2.175)  and  (2.176),  we  find 


u;r 


(2.1  S3) 
(2.184) 

These  equations  of  motion  can  be  written  in  the  frequency  domain  by  defining 


r^  +  —r,-r^;r,  =  n,. 
Vr 

6,  +  —0,  +  —  O,. 

Qe 


r,(t) 


r.  u;  e 


(2.185) 
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a.(0  =  a.(u;)e'“", 

e,(t)  =  e,{u)e'^\ 
a.(/)  =  a.(a;)e“"'. 


(2.186) 

(2.187) 

(2.188) 


Sub«*itutinp  thpsf*  deflnitions  into  our  equations  of  nioliun,  and  solving  for  the  ratios, 
r,(a;)/a,(c4;)  and  ^,(tj)/o,(a;),  we  find 


If  ax  — 

fiay  = 
Haz  = 

lice.  = 
iioe,  ^ 
Has.  = 


Tt{^) 

1 

+  iu>fU}/Q.r  —  U)'^ 

r„(u;) 

1 

ay{u) 

LJ^  -f-  iiJrUjfQr  -  ijp- 

T^{u) 

1 

0^(02) 

T  ii^r'^IQr  -  LtP-  ’ 

e,(u) 

1 

Ox(t^) 

Jl  +  iLJeLjjQg  - 

1 

Oj,(u;) 

+  iugulQe  -  u>'^ 

1 

Ox(w) 

uj  +  kOgijfQg  - 

(2.189) 

(2.190) 

(2.191) 

(2.192) 

(2.193) 

(2.194) 


The  transfer  functions  relating  acceleration  to  displacement,  Eqs.  (2.189)  through  (2.194). 
and  displacement  to  output  current,  Eqs.  (2.160)  through  (2.165),  can  be  combined  to 
define  the  transfer  functions  relating  acceleration  to  output  current: 

ft  _  , _  ,  i^S  I  ,r,  lnx^ 

**ax\  —  /  X  ~  IfCOSiJrx^  j  I  c.  T  2  i  '  in  2'  (2.196) 

rtx(^)  L>SQ  +  6Ls  ^ 

Hff,,  =  ^  =  IrCOSUr,t--{~ -  -  \  - 2’  (2196) 

LsQ  +  oLs  +  lUrU>IQf  - 

Tf  _  r  ■  ^5  1  ^r,■7\ 


ir  —  ^xv~;  T  ,  “j 

..  Oi«  =  ,  \  ~  /r<^05u7rj.t  i  CJ  2  i  '  //n  2’ 

ai(u)  LsQ  +  6L5  Uf  +  tUJri^lQr  -  U!^ 

U  -  _  r _  .  1 

LsQ  +  oLs  +  lUrU>IQf  - 

IJ  -  *^(‘^)  _  r .  ,  1 

iia.i  —  /  \  ~  ifCOSijJr^l  ,  n  r  2  1  ’  //O  2' 

a,[ijj)  i'SQ  T  u).  +  liJrijJ I Qf  —  ui^ 

u  _  ,  ,  fA5  1 

-  0,(0;)  "  +  6L5  +  r^e^lQo  -  0^2  ' 

fj  _  *ev(‘^)  _  r  .  1 

-  o,(u;)  -  ^  6Z.S  u;2  +  ,u;«u;/g,  -  u;^  ' 

ft  _  r  ,  <^A5  1 

-  a,(u>)  + 


(2.195) 

(2.196) 

(2.197) 

(2.198) 

(2.199) 

(2.200) 


2.5  Fundamental  Noise  of  the  Accelerometer 


The  fundamental  noise  of  the  SSA  can  be  described  in  terms  of  the  equivalent 
acceleration  noise,  and  ,  applied  to  the  proof  mass.  'I'his  equivalent  input  accelera- 
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tion  noise  is  generated  by  two  separate  processes  in  the  SSA.  The  first,  Brownian  motion 
noise,  actually  generates  acceleration  noise  in  the  form  of  thermodynamic  fluctuations. 
The  second,  SQUID  amplifier  noise,  adds  noise  to  the  measured  signal,  which  is  reflected 
in  the  output  as  an  equivalent  input  acceleration. 

Since  the  processes  in  the  SSA  are  stationary,  we  can  equate  the  auto  corrchuion 

+  r))  (2.201) 


of  a  signal  from  the  SSA,  x(<),  to  a  time  average: 

RixiT)  =  lini  ^  f  x{t)x'‘{t  T)dt.  (2.202) 

T— oo  21  J-T 

The  Wiener-Khintchine  theorem  [14,15]  allows  us,  in  turn,  to  relate  this  to  the  power 
spectrum  density: 

=  r  RxT{r)e-'^"dT.  (2.203) 

J  «-oo 

Rxxir)  ^  ^  r  Sx{oj)e'^^cL;.  (2.204) 

2ir  y_oo 

2.5.1  Transducer  Brownian  Motion  Noise 


The  Brownian  motion  of  the  SSA  at  a  temperature  T  appears  at  the  proof 
mass  as  three  linear  and  three  angular  noise  terms,  aj(u>)  and  of(w),  respectively.  By 
applying  aj{u)  and  aj{u)  to  the  transfer  functions  relating  the  position  and  orientation 
of  the  SSA  to  acceleration,  Eqs.  (2.189)  through  (2.194),  we  can  relate  the  position  and 
orientation  of  the  SSA  to  the  Brownian  motion  noise.  This  will  allow  us  to  compute 
the  average  potential  energy  of  each  mode,  and  in  turn,  by  applying  the  equipartition 
theorem,  compute  the  spectral  density  of  the  Brownian  noise. 

The  transfer  functions  relating  the  position  of  the  SSA  to  the  applied  Brownian 
motion  noise  are,  from  Eqs.  (2.189)  through  (2.194). 


r,{u;) 

d,(u>) 


aji^) 

+  iulijJr  j Qr  ~ 

UJg  T  iu)u)gfQg  —  u/'*' 


(2.205) 

(2.206) 


The  thermal  energy  in  each  of  the  linear  modes  of  the  SS.4  can  now  be  computed: 


ksT  =  om;2(r?(<)) 


4^  f  Sr,(^)d^ 

J  — oo 

■^tt—  /  - — ^ 

^TT  y  — oc>  Icy?  “I"  lOJLJ 


'In  y_oc  |cy?  +  iuLJr/Qr  -  cy^j 

where  kg  is  Boltzmann’s  constant  and  we  have  used  the  fart,  fr^rn  Kqs.  (2.205),  (2.202), 
and  (2.203),  that 

^  5'aT(cy) 

5,  (cy)  = - = - r-  (2.20S) 

icy?  +  ;cyu- /r;,  -cv-^!' 

Since  the  thermal  fluctuations  are  white  in  character,  i.s  a  constant,  and  we  can 

I 

perform  the  integral  immediately.  We  find 


,  mQr  „ 


(2.200) 


,  2kBTu;, 


(2.210) 


This  equation  is  the  mechanical  analog  of  the  Nyquist  e(iuation.  I'lie  power  si)ectra! 
density  for  each  of  the  angular  modes  of  the  SSA  can  similarly  be  computed; 


5,,r  = 


'IksTug 


(2.211) 


2.5.2  SQUID  Amplifier  Noise 

A  noise  model  of  the  rf  SQUID  can  be  found  in  [6.7].  In  this  model,  the 
optimum  source  impedance  is  much  smaller  than  the  input  impedance  ujLsq-  Because 
of  this  impedance  mismatch,  the  noise  of  the  SQUID  is  dominated  by  current  noise, 
5/(cy),  and  can  b«  characterized  by  an  “input  energy  sensitivity": 

Es{i^)  =  ~LsQS!(^).  (2.212) 

2.5.3  Potential  Sensitivity 

The  output  of  the  SSA  is  modulated  by  the  applied  carrim  frequency  in  each  of 
the  six  degrees  of  freedom  and  each  of  these  signals  must  be  separately  demodulated  by  a 
lock-in  amplifier.  This  process  is  represented  schematically  in  Fig.  2.9,  for  the  Tj.  degree 
of  freedom.  In  order  to  fully  understand  how  each  fundamental  noise  source  contributes 
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Figure  2.9.  A  schematic  illustrating  detection  of  acceleration  in  the  SS.4: 
Tx  degree  of  freedom. 


to  the  total  noise  of  the  SSA,  we  will  travel  through  each  stage  of  the  process  for  the 
degree  of  freedom. 

Any  acceleration  signal,  drives  the  the  SSA  must  compete  with  the 

Brownian  motion,  a'^(t).  Thus  the  combined  power  spectral  density  which  appears  at 
the  input  to  the  block  diagram  in  Fig.  2.9  is 


5aju;)  +  SaT. 


(2.2i:3) 


This  spectral  density  is  translated  by  the  frequency  response  of  the  SSA  into  a  displace¬ 
ment  spectral  density: 


SrAu)  = 


SaA>^)  + 


(2.214) 


|cj^  -f  iuUrlQr  - 

This  signal  is  modulated  and  appears  at  the  input  to  the  SQUID  amplifier.  Modulation 
of  the  Tx  signal,  Eq.  (2.160),  relates  the  spectral  density,  S,^.  at  the  input  to  the  SQUID, 
to  the  spectral  density  Sr^'. 


SiA^^) 


1 


IrA 


ri\S 


4  [LsQ-hGLsJ 
Combining  this  with  Eq.  (2.214).  we  find 


[.Sr_j(a-'  —  )  +  •Srj(^'  +  )]  ■ 


(2.215) 


5’.x(‘^)  =  7 


/rA< 


4  \  Lhq  -b  67/5 


k?  +  '(w’  -  -v,  ,i-v/Qr  --  (w  -  )‘| 


+ 


+  U-Vi  )  +  .5’g7 


+  5’/, 


(2.216) 


Mr  T  y^'rlQr  —  (al  +  )^|^_ 

where  we  have  added  the  input  noise  spectral  den.sity  from  the  SQUID  amplifier,  .5/. 

4’his  signal  is  then  amplified  by  the  SQUID  and  (h'modulated.  The  power 
spectral  density,  just  before  the  low  pass  filter  in  Fig.  2.9.  is 

‘Si(^’)  =  —  wVi  )  +  +  w'r,  )]  . 


(2.217) 
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where  G  is  the  gain  of  the  SQUID  in  V/A.  Substituting  Eq.  (2.216),  we  find 


1  /  Ir^s  j  )  +  5gr 

4  yi-sg  +  67/5/  |a;2  +  i(ai  -  2ur^)ujJQr  -  {u>  -  2t^r 

+  2Wr,)  +  +  ■S'aJ'  11 


S.A^)  =  -GU2S,  +  - 


_ t- - ;*£ - -  -I-  4 - ^-2 — : - -  \ . 

\ul  +  l{u  +  ‘hjSrA^rlQr  -  (tc>  +  2LJrAA^  \‘-^t  +  ^'^rlQr  “  j 

(2.218) 

Note  that  a  4,  not  a  2,  appears  in  front  of  the  last  term  inside  the  brackets,  because 
the  signals  are  adding  coherently.  Since  the  resonances  of  the  SSA  are  of  a  much  lower 
frequency  than  0;^^,  the  first  two  terms  within  the  brackets  on  the  right  hand  side  of  Eq. 
(2.218)  can  be  neglected,  and  the  final  output  spectral  density  is 


■5u.(a;)  =  <^25/  + 


/rA5  r  +  1 

^sg  +  6^5/  +  iuIWr/Qr  -  J 


(2.219) 


If  we  represent  the  SQUID  amplifier  noise  as  an  equivalent  acceleration  noise, 
af*^,  acting  on  the  proof  mass,  then  Eq.  (2.219)  can  be  rewritten: 


4  [LsQ  +  Us)  [  \u}  +  iwUrlQr-^^?  V 


(2.220) 


where 


5  5q(w)  =  2 


_  /i^+6lsVi  ,2 


+  iuJUJrlQr  -  Si- 


(2.221) 


Thus  the  minimum  detectable  acceleration  spectral  density  is.  from  Eq.  (2.220), 


>  S  Sq{u})  +  S^T 


^  ,  (  LsQ  +  QLs\^  \  2  ,  ■  lr^ 

~  4  ^ - )  l’’  ^  ^^^r/Qr  -  I 


2  Es  2kBTu>r 
ksQ  mQr 


,  (2.222) 


where  we  have  substituted  Eqs.  (2.221),  (2.210),  and  (2.212). 

The  power  spectral  density  extends  over  both  positive  and  negative  frequencies, 
and  in  order  to  compute  the  noise  power,  Pq,  in  any  frequency  band,  D,  we  must  include 
both  the  negative  and  positive  frequencies.  Thus,  it  is  convenient  to  define  a  one  sided 
power  spectral  density,  P{>^)- 


P(u>)  =  [.8(u;)  +  S(  -a;)] , 


and  the  power  in  a  frequency  band,  Q,  is 


Pu=  I  Pi2iTf)df. 
Ju 


(2.223) 


(2.224) 
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The  minimum  detectable  acceleration  spectral  density,  Sa^(uj),  is  then  equivalent  to 


PaA^)  =  8  - 


2\'^  Es  AksTur 


V  /.As  ;  I  ^  I  /^SQ  mQr 

For  frequencies  much  less  than  the  resonances  of  the  SSA,  this  simplifies  to 


D  _  o  /" ^SQ  +  6/,s\^  4  Es  ,  AksTur 

--H~7;a7  ) 


(2.225) 


(2.226) 


This  can  be  written  in  terms  of  an  energy  coupling  coefficient,  Pr- 


\mu^rl 


(2.227) 


where  (  )t  denotes  a  time  average  over  one  period  of  the  carrier  frequency  Ur^.  Substi¬ 
tuting  Eq.  (2.160)  into  our  expression  for  /?r,  we  find 


3  =-(  1  LAR. 

2  \Lsq  +  6Ls  j  mu)'} 


(2.228) 


Substituting  this  into  the  expression  for  the  minimum  detectable  acceleration  spectral 
density,  Eq.  (2.226),  gives 


p  j. 

mKQr  dr 


(2.229) 


Comparing  this  to  the  minimum  expression  for  the  noise  in  the  gravity  gradiometer,  it 
is  seen  that  the  noise  from  the  SQUID  contributes  twice.  This  is  a  consequence  of  the 
demodulation  which  folds  the  two  separate  sidebands  of  the  modulated  signal  upon  each 
other,  so  that  the  noise  in  each  sideband  is  added  together.  If  we  had  used  six  separate 
SQUIDs  in  a  more  traditioncil  DC  detection  scheme,  we  would  not  have  suffered  this 
penalty,  but  we  would  have  lost  all  the  benefits  gained  from  modulation. 

The  coupling  coefficient,  dr  can  be  rewritten  in  a  more  compact  form.  If  we 
remove  the  direct  dependence  on  /^  and  u>r,  by  substituting  Eqs.  (2.155)  and  (2.172)  into 
Eq.  (2.228),  we  find 


Lsq  +  S/'S  yki  +  ks / 

Thus  we  see  that  dr  can  vary  between  0  and  1/(1  +6/,s//'Sg)-  If  Brownian  noise  in  the 
SSA  is  dominated  by  SQUID  noise,  is  minimized  when  ks  =  k^,  or 


1  Lsq 

2  /s£j  +  6/s 


(2.231) 
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This  optimal  coupling  could  be  improved  by  inserting  a  superconducting  transformer 
between  the  bridge  circuitry  and  the  SQUID.  This  would  match  Lsq  +  6Ls  to  Lsq,  and 
/?r  optimal  would  be  1/2. 

A  similar  analysis  of  the  two  other  linear  degrees  of  freedom  gives 


m  \  Qr  '  0r 
Likewise,  analysis  of  the  three  angular  degrees  of  freedom  gives 


where 


.  /  LsqijX  ^  1 

*  \  /,  2 


(2.232) 


e  /ksT  u;eEs\ 

\Qe  09  )' 

(2.233) 

[gcAs  \  Lsq 

Lsq  +  6Ls  J  luj 

(2.234) 
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Chapter  3 

Construction  of  the 
Superconducting  Six- Axis 
Accelerometer 


In  the  previous  chapter  we  discussed  the  theory  of  the  SSA.  In  this  chapter  we 
will  describe  the  actual  realization  of  this  theory  in  a  concrete  design.  In  the  first  section, 
we  will  describe  the  main  mechanical  components  that  comprise  the  SSA:  the  proof  mass, 
coil  forms,  coil  form  holders,  precision  mounting  cube,  transformers,  and  superconducting 
coils.  In  the  second  section,  we  will  describe  the  superconducting  circuitry.  In  the  third 
section,  we  will  examine  the  design  of  the  cryostat  and  the  suspension  hardware  that 
supports  the  SS.A  in  the  cryostat.  In  the  fourth  and  final  section,  we  describe  the 
computer  system,  data  acquisition,  and  interface  electronics  used  to  control  the  SSA. 


3.1  Accelerometer  Hardware 

3.1.1  Mechanical  Components 

The  primary  hardware  that  makes  up  the  SS.\,  consists  of  four  main  compo 
nents:  the  proof  mass,  the  titanium  coil  forms  which  hold  the  superconducting  coils  in 
close  proximity  to  the  proof  mass,  the  coil  form  holders  which  align  the  coil  forms,  and 
a  precision  mounting  cube  which  holds  the  coil  form  holders.  An  exploded  view  of  the 
SSA  is  shown  in  Fig.  3.1. 
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The  proof  mass  for  the  SSA  is  constructed  from  pure  niobium  (Nb).  It  consists 
of  four  separate  plates  which  fit  tightly  together  to  form  the  final  proof  mass.  An 
exploded  view  of  this  is  shown  in  Fig.  3.2.  The  machining  of  the  proof  mass  proceeded 
in  four  steps.  In  the  first  stage  the  four  pieces  were  cut  from  a  Nb  plate.  Each  piece 
was  left  slightly  oversized  in  order  to  maintain  the  precision  of  the  finaJ  proof  mass.  The 
pieces  were  then  heat  treated  in  a  high  vacuum  oven  at  1800°C  for  4  hours.  Next  all 
dimensions  except  the  slots  were  finished  using  a  computer  controlled  milling  machine. 
The  slots  were  then  finished  using  an  electron  discharge  machine  (EDM).  The  pieces 
were  then  polished  flat  to  a  final  thickness  of  2.39  mm  and  assembled  into  the  final  proof 
mass.  A  small  amount  of  GE  7031  varnish  (General  Electric  Co.,  Schenectady,  New 
York)  was  used  to  lock  the  pieces  in  place.  The  final  dimensions  of  the  proof  mass  are 
2.784  cm  on  a  side.  The  proof  mass  is  orthogonal  and  flat  to  a  part  in  10^. 

The  proof  mass  is  surrounded  by  eight  titanium  alloy  cubes  which  hold  the 
superconducting  coils  in  close  proximity  to  the  proof  mass.  Each  of  the  eight  cubes  has 
three  surfaces  which  face  the  proof  mass,  and  three  surfaces  which  do  not.  Each  of  the 
surfaces  which  face  away  from  the  proof  mass  has  a  recessed  pocket  in  which  small  strain 
relief  clamps  are  mounted.  These  clamps  protect  the  superconducting  wires  that  connect 
to  the  superconducting  coils.  Each  of  the  surfaces  which  face  the  proof  mass  has  two 
flat  circular  regions  which  contain  one  levitation  coil  and  one  sensing  coil,  respectively. 
An  isometric  drawing  of  a  coil  form  is  shown  in  Fig.  3.d.  Titanium  was  chosen  as  the 
material  for  the  coil  forms  because  of  its  strength  and  machinability.  It  was  also  felt 
that  a  conducting  coil  form  would  decrease  the  susceptibility  of  the  superconducting 
circuitry  to  rf  noise,  since  eddy  current  losses  add  an  additional  loss  mechanism  at  high 
frequencies. 

The  coil  forms  are  held  in  place  by  a  set  of  six  cylindrical  Nb  coil  form  holders 
which  mount  inside  a  precision  10.16  cm  titanium  cube.  These  six  coil  form  holders 
mesh  together  inside  the  titanium  cube  in  such  a  way  as  to  form  a  smaller  hollow  Nb 
cube  inside  the  titanium  cube.  The  eight  coil  forms  are  mounted  on  the  inside  of  this 
hollow  Nb  cube.  Each  coil  form  is  attached  to  three  of  the  six  Nb  coil  form  holders. 
This  is  possible  since  the  six  coil  form  holders  form  the  six  sides  of  the  hollow  Nb  cube, 
and  each  coil  form  rests  in  one  of  the  eight  corners.  Each  of  the  recessed  pockets  in  the 
coil  forms  is  matched  by  a  corresponding  hole  in  the  coil  form  holder.  Wiring  from  the 
superconducting  coils  on  the  coil  forms  is  routed  through  these  holes  to  the  large  circular 


Figure  3.4.  An  isometric  view  of  the  six  coil  form  holders. 


^.0''‘^cts  i^i  rhe  coil  form  holders.  All  the  superconducting  joints  in  the  SSA  are  mounted 
in  these  large  circular  pockets.  Small  holes  allow  the  superconducting  wires  to  be  routed 
from  one  pocket  to  another. 

The  cylindrical  pocket  at  the  center  of  each  coil  form  holder  is  used  to  house  a 
feedback  transformer.  An  isometric  view  of  the  six  Nb  coil  form  holders,  as  they  would 
appear  inside  the  precision  mounting  cube,  is  shown  in  Fig.  3.4. 

The  precision  titanium  cube  which  houses  the  SS.\  is  shown  in  Fig.  3. -5.  Its 
surfaces  are  flat  and  mutually  orthogonal  to  better  than  2  parts  in  10'’.  Each  of  the  six 
surfaces  on  the  mounting  cube  contains  a  tapped  bolt  circle  (not  shown)  which  match 
the  clearance  holes  in  the  accelerometers  wiiich  make  up  the  superconducting  gravity 
gradiometer.  Presently  these  bolt  circles  are  used  to  mount  .Vb  covers  over  the  pockets 
of  circuitry  in  the  SSA.  but  when  the  SSA  is  integrated  with  the  gravity  gradiometer, 
these  covers  will  be  replaced  by  the  six  accelerometers  that  make  up  the  superconducting 
gravity  gradiometer. 


Figure  3.5.  An  isometric  view  of  the  precision  mounting  cube. 


Each  corner  of  the  mounting  cube  has  been  cut  perpendicular  to  a  line  passing 
through  opposite  vertices  of  the  cube.  One  of  these  corners  is  used  to  mount  the  SSA 
in  the  cryostat  on  the  suspension  structure.  Current  and  control  leads  enter  the  SSA 
through  the  other  seven  corners. 

3.1.2  Superconducting  Coils 

Each  surface  of  the  coil  form  facing  the  proof  mass  contains  two  concentric 
superconducting  coils.  The  outer  coil  is  used  to  provide  levitation  and  feedback  to 
the  proof  mass.  The  inner  coil  is  used  to  sen.se  the  position  of  the  proof  mass.  Each 
coil  is  a  spiral  pancake  coil  wound  from  76  /rm  diameter  niobium-titanium  (Nb-Ti) 
superconducting  wire.  Originally  we  constructed  the  coils  in  the  SSA  using  pure  Nb 
wire,  but  recent  changes  in  the  manufacturing  process  for  pure  Nb  wire  created  reliability 
problems  which  we  art  irrently  attempting  to  understand.  Since  the  SSA  is  constructed 
from  48  superconducting  coils,  even  a  small  failure  rate  is  unacceptable  and  we  switched 
to  tougher  Nb-Ti  wire.  A  small  amount  of  flux  creep  has  been  observed  in  Nb  Ti  coils 
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[9]  and  this  introduces  additional  noise  into  any  superconducting  circuitry  containing 
Nb-Ti  coils.  This  is  not  yet  a  factor  in  the  SSA  since  our  noise  is  still  dominated  by 
seismic  disturbances  in  the  laboratory. 

The  Nb-Ti  wire  is  insulated  with  ML  insulation  and  has  an  overall  diameter 
of  100  fim.  Before  the  Nb-Ti  wire  is  wound  on  the  titanium  coil  form,  the  surface  is 
carefully  sanded  and  cleaned  of  all  dirt  and  oil.  The  outer  most  coil  is  wound  by  placing 
a  clear  highly  polished  Plexiglas  backing  plate  on  top  of  the  coil  form.  This  backing 
plate  rests  on  top  of  the  surface  of  the  inner  coil  and  is  held  in  place  with  a  boll  during 
the  winding  process.  There  exists  a  small  (200  /xm)  step  between  the  outer  and  inner 
coil.  This  small  step  prevents  the  backing  plate  from  making  contact  with  the  outer 
coil  winding  su-face,  and  leaves  room  for  a  two  layer  spiral  coil.  A  thin  layer  of  vacuum 
grease  is  applied  to  the  surface  of  the  inner  coil  to  prevent  epoxy  from  bonding  to  the 
surface  of  the  inner  coil  during  the  winding  of  the  outer  coil.  A  small  radial  groove  from 
the  center  of  the  coil  form  to  the  outside  allows  the  insertion  of  two  small  hollow  Teflon 
tubes.  The  first  tube  runs  from  the  center  of  the  inner  coil  outwards  and  allows  us  to 
wind  the  inner  coil  after  the  outer  coil  is  finished.  The  second  tube  guides  the  start  of 
the  Nb  Ti  wire  used  in  the  outer  coil.  Winding  starts  from  tfip  center  of  the  outer  roil 
and  ends  after  a  counted  number  of  turns.  A  transparent,  low  viscosity,  epoxy,  TRA- 
CAST  BB3002  (TRA-CON,  Inc..  Medford,  Massachusetts)  is  used  to  bond  the  Nb-Ti 
wire  to  the  coil  form.  .A.fter  the  epoxy  is  set,  the  center  bolt  is  removed  and  the  backing 
plate  is  machined  off,  leaving  a  flat  spiral  coil  consisting  of  two  layers.  The  coil  winding 
procedure  used  here  is  a  modification  of  that  aescribed  in  [10, 3j. 

After  the  outer  coil  is  finished  and  inspected,  the  inner  coil  is  wound.  A  hollow 
titanium  pin  with  a  small  step  (100  /xm)  is  inserted  into  the  center  hole  of  the  coil  form. 
The  small  step  prevents  epoxy  from  running  down  the  central  hole.  A  small  amount  of 
vacuum  grease  is  applied  to  the  pin  to  prevent  a  bond  from  forming  between  the  pin  and 
the  epoxy.  A  second  clear  Plexiglas  backing  plate  with  a  center  hole  of  diameter  equal 
to  the  titanium  pin  is  mounted  on  the  pin.  Careful  adjustment  of  the  spacing  between 
the  backing  plate  and  the  surface  of  the  coil  form  matches  this  gap  to  the  diameter  of 
the  Nb-Ti  wire.  A  small  bolt  inserted  through  the  hollow-  titanium  pin  is  used  to  hold 
this  spacing  constant.  Winding  proceeds  from  the  center  of  the  inner  coil  and  ends  after 
a  counted  number  of  turns.  After  the  epoxy  is  set,  the  center  pin  is  removed  and  the 
backing  plate  is  machined  off,  leaving  a  flat  spiral  coil  consisting  of  a  single  layer. 
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Each  of  the  pairs  of  Nb-Ti  wire  from  the  two  coils  is  then  inserted  in  color  coded 
Teflon  tubes.  Each  Teflon  tube  is  anchored  to  a  small  stress  relieving  clamp  mounted  on 
the  adjoining  surface  of  the  coil  form  that  faces  away  from  the  proof  mass  surface. 

After  all  parts  are  cleaned  and  inspected,  the  complete  SSA  is  assembled.  A 
photograph  of  the  assembly  is  shown  in  Fig.  3.6.  The  large  number  of  leads  from  the 
superconducting  coils  can  be  clearly  seen. 


3.2  Superconducting  Circuitry 

The  superconducting  joints  used  in  the  SSA  are  formed  by  spot  welding  two 
or  more  Nb-Ti  wires  to  Nb-Ti  foil.  Each  small  section  of  foil  has  two  or  more  small 
sections  of  Teflon  tubing  mounted  at  one  end.  These  tubes  keep  the  wire  from  bending 
at  the  spot  weld  and  provide  strain  relief.  Before  the  spot  weld  is  formed,  each  wire 
is  crushed  flat  to  approximately  half  its  original  diameter.  This  increases  the  surface 
area  of  the  wire  in  contact  with  the  foil.  Both  the  Nb-Ti  wires  and  the  Nb-Ti  foil  are 
polished  in  a  diluted  chemical  solution  of  equal  parts  85%  H3PO4,  70%  HNO3,  and  49% 
HE,  and  cleaned  with  distilled  water  to  remove  any  residual  acid.  The  upper  and  lower 
electrodes  used  in  the  spot  welder  are  made  of  beryUium  copper.  The  upper  electrode 
makes  contact  with  the  flattened  wire  and  is  rectangular  in  shape,  approximately  1  mm 
by  5  mm.  The  lower  electrode  is  cylindrical  in  shape,  approximately  5  mm  in  diameter. 
The  contact  pressure  and  energy  settings  of  the  spot  welder  are  adjusted  so  that  a  test 
specimen,  when  viewed  under  a  microscope,  shows  a  clear  fusion  of  the  wire  and  the  foil 
with  little  or  no  evaporation  of  material.  Superconducting  joints  formed  in  this  manner 
consistently  have  critical  currents  in  excess  of  8  A  through  many  thermal  cycles. 

The  superconducting  feedback  transformers  used  in  the  SSA  are  cylindrical  in 
shape,  with  a  75  /im  thick  layer  of  brass  separating  the  primary  and  the  secondary.  This 
thin  brass  layer  forms  a  barrier  to  high  frequency  interference,  and  allows  a  direct  con¬ 
nection  of  the  feedback  electronics  to  the  SSA.  Six  feedback  transformers  are  employed 
in  all.  These  are  mounted  inside  the  cylindrical  cavities  in  the  six  Nb  coil  form  holders. 

Six  additional  tuned  transformers  are  used  to  connect  the  six  inductance  bridges 
to  their  room  temperature  oscillators.  Each  transformer  steps  up  the  current  supplied 
by  the  room  temperature  oscillator  by  a  factor  of  six.  These  transformers  are  cylindrical 
in  shape  and  enclosed  in  2.54  cm  diameter  Nb  tubing.  A  thin  brass  layer  separates  the 
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primary  and  the  secondary,  and  forms  a  barrier  to  high  frequency  interference.  This  is 
extremely  important  since  the  secondary  of  each  transformer  is  attached  directly  to  the 
sensing  circuitry,  and  very  low  levels  of  rf  interference  will  cause  the  SQUID  control  loop 
to  unlock.  Each  transformer  has  a  polystyrene  capacitor  attached  across  the  primary. 
This  allows  the  transformer  to  resonate  at  a  particular  frequency,  forming  an  additional 
noise  filter  for  frequencies  other  than  the  resonance  frequency,  and  providing  additional 
current  gain  at  that  frequency.  The  effectiveness  of  this  scheme  is  demonstrated  by  the 
ability  of  the  SQUID  to  stay  locked.  Each  of  these  transformers  is  mounted  on  the 
underside  of  the  top  of  the  vacuum  can  which  houses  the  SSA. 

Shielded  heat  switches  are  made  from  tiny  1  /8  W,  500  Cl  carbon  resistors. 
These  are  specially  selected  carbon  resistors  whose  resistance  increases  dramatically  as 
they  approach  liquid  helium  temperatures.  This  provides  passive  feedback  when  the 
heat  switch  is  driven  by  a  current  source.  Each  resistor  is  insulated  and  surrounded 
by  a  very  thin  layer  of  lead  foil.  This  is  then  soldered  to  a  lead  tubing  which  contains 
the  leads  that  attach  to  the  resistor.  This  shielding  prevents  rf  frcqucncj  f  cm  entering 
the  c’lperconducting  circuitry  through  the  leads  connected  to  the  heat  switch.  A  Nb-Ti 
wire  is  then  folded  in  half  and  wound  noninductively  around  the  lead  shielded  resistor. 
This  wire  is  held  in  place  with  a  small  amount  of  GE  7031  varnish.  The  completed  heat 
switch  is  attached  to  a  small  copper  block  that  contains  a  round  copper  post  that  is 
tinned  with  solder.  The  Nb-Ti  wire  is  then  wound  several  times  around  this  post  and 
held  in  place  with  additional  GE  7031  varnish.  The  copper  post  serves  as  a  heat  sink 
and  insures  that  the  only  portion  of  the  superconducting  circuitry  that  is  warmed  above 
the  transition  temperature  is  the  Nb-Ti  wire  directly  in  contact  with  the  heat  switch. 
These  shielded  heat  switches  are  then  mounted  inside  the  large  circular  pockets  in  the 
Nb  coil  form  holders.  In  all,  sixteen  heat  switches  are  used  in  the  SS.\. 

3.2.1  Levitation  Circuitry 

The  levitation  circuits  used  in  the  SSA  to  support  and  control  the  x  linear  and 
angular  degrees  of  freedom  are  shown  in  Fig.  3.7,  wliere  all  heat  switches,  superconduct¬ 
ing  joints,  and  current  leads  are  shown  explicitly.  The  y  and  ^  degrees  of  freedom  have 
identical  (in  form)  levitation  circuits.  Flach  of  the  superconducting  Joints  is  represented 
by  a  black  dot.  Each  heat  switch  is  shown  schematically  as  a  resistor  in  close  proximity 
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to  a  superconducting  lead.  Each  heat  switch  has  associated  with  it  a  pair  of  current 
leads  which  are  used  to  store  persistent  currents  in  the  circuitry. 

All  current  leads  in  the  SSA  are  constructed  with  125  /rm  Cu  clad  Nb--Ti 
wire.  The  copper  cladding  reduces  the  risk  of  driving  the  current  supply  leads  normal, 
and  provides  additional  heat  capacity  that  protects  the  leads.  A  programmable  curnuit 
supply,  controllable  though  the  GPIB  interface,  is  used  to  supply  external  current  to  these 
current  leads.  If  any  portion  of  the  superconducting  circuit  is  driven  normal  while  storimT 
currents,  the  external  current  supply  would  damage  the  superconducting  circuitry.  In 
order  to  prevent  this,  an  automatic  protection  circuit  monitors  the  voltage  across  the 
current  supply  leads.  Any  voltage  signal  which  exceeds  a  preset  voltage  threshold  triggers 
a  current  shunt. 

The  r~  levitation- feed  back  circuit  is  less  complex  that  the  levitation-feedback 
circuit.  Only  a  single  heat  switch,  hsl.  is  u.sed  lo  disrupt  tlie  superconducting  loop.  The 
storage  of  a  persistent  current  in  the  superconducting  loop  proceeds  in  several  steps. 
Initially  there  is  no  current  stored  in  the  loop.  The  external  current  leads,  /],  are  then 
driven  with  a  steady  current  of  0.5  A.  At  this  point,  all  the  current  from  the  external 
current  supply  driving  I\  is  passing  through  the  heat  switch,  hsl  is  then  pulsed  on  for 
approximately  40  msec.  This  drives  the  short  path  between  the  two  superconducting 
joints  normal  and  forces  the  current  to  flow  through  the  four  inductors  which  levitate 
the  proof  mass  in  the  direction.  Approximately  10  msec  after  the  current  driving 
the  heat  switch  is  removed,  the  heat  switch  becomes  suj>erronducling  once  again.  This 
traps  the  current  in  the  r^.  levitation- feed  back  loop.  This  process  is  repeated  at  higher 
and  higher  current  settings  until  the  desired  amount  of  flux  is  trapped  in  the  circuit. 
The  bridge  is  carefully  balanced  in  an  iterative  process  which  is  described  in  the  next 
chapter. 

The  dr  levitation-feedback  circuit  contains  two  heat  switches,  hs2  and  hs3.  and 
two  pairs  of  current  supply  leads,  I-i  and  Initially,  current  i>  stored  about  the  outci 
loop  of  the  circuit.  This  is  done  by  i)ulsing  hs2  ami  h.sd  on  while  su[)plying  current  to 
This  traps  flux  in  the  circuit  so  that  —  f2(\-  t  his  provides  an  even  downward 

pressure  on  the  proof  mass  in  the  — .r  direction.  If  a  currmil  is  now  supplied  through 
and  hs3  is  pulsed  on,  the  current  tra[>|)ed  in  each  side  of  the  circuit  will  be  such  that 
—  h  —  and  =  I2  +  h/'*  4  his  imbalance  in  current  apjilies  a  torrpie  atmut 

the  axis  of  the  proof  rmuss.  In  practice,  /(  is  used  to  a[)i)ly  a  small  tonjue  to  the  proof 
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mass  which  balances  the  sensing  bridge. 

3.2.2  Sensing  Circuitry 

The  sensing  circuit  used  in  the  S3A  is  shown  in  Fig.  3.8,  where  all  heat  switches, 
superconducting  joints,  and  current  leads  are  shown  explicitly. 

The  six  heat  switches,  hs4-hs9,  are  used  to  discharge  any  flux  that  may  be 
trapped  in  the  six  sensing  bridges.  The  seventh  heat  switch,  hslO,  is  used  to  discharg«' 
any  flux  that  may  be  trapped  in  the  superconducting  loop  containing  the  SQUID  input 
coil.  This  heat  switch  also  serves  to  protect  the  SQUID  from  any  input  current  spikes 
that  may  be  induced  in  the  sensing  circuitry  while  currents  are  being  stored  in  the 
levitation  circuitry. 

All  leads  from  the  sensing  inductances  are  of  equal  length,  and  both  sides  of 
each  sensing  bridge  are  identiccil  up  to  the  small  inductance  added  by  each  heat  switch. 
This  symmetry  is  important  as  the  additional  leads  in  the  sensing  circuitry  generate 
stray  inductance  which  can  imbalance  the  sensing  bridge.  Any  imbalance  forces  us  to 
displace  the  proof  mass  away  from  equilibrium  in  order  to  rebalance  the  sensing  bridge. 
Presently  the  heat  switches  in  the  SSA  sensing  circuit  only  lie  on  one  side  of  each  bridge. 
It  is  possible  to  construct  an  eight  wire  heat  switch  which  would  present  a  symmetrical 
set  of  stray  inductances  to  the  sensing  bridge,  but  the  present  choice  was  made  for  its 
simplicity. 

The  six  current  leaxls,  arc  driven  by  six  room  temperature  oscillators.  .A. 

polystyrene  capacitor  is  connected  across  the  primary  of  each  transformer  to  form  a  tank 
circuit.  The  tank  circuit  strongly  attenuates  noise  above  its  resonance  frequency  and. 
when  driven  on  re.sonance,  increases  the  driving  current  through  the  inductance  bridge. 
The  capacitors  are  selected  so  that  each  tank  circuit  resonates  at  a  separate  frequency 
with  the  lowest  at  a  frequency  of  500  IIz,  and  the  highest  at  a  frequency  of  2500  IIz. 
The  resonant  frequencies  lie  equally  spaced  throughout  this  band. 

The  output  of  the  SQUID  amplifier  (Model  2000,  Quantum  Design,  Inc.,  San 
Diego,  CA)  is  high-pass  filtered,  and  fed  into  six  separate  lock-in  amplifiers  (ModpI 
3961B,  Ithaco,  Inc.,  Ithaca,  NY).  Each  lock-in  amplifier  is  tuned  to  one  of  the  six  carrier 
frequencies  and  one  of  the  six  .separately  modulated  signals  is  demodulated  and  recovered 
at  the  output  of  the  lock-in. 
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3.3  Construction  of  the  Experimental  Test  Apparatus 

In  this  section,  we  describe  the  experimental  test  apparatus  which  houses  the 
SSA:  the  thermal  and  mechanical  design  of  the  cryostat,  the  electrical  isolation  of  the 
control  leads,  th«>  vibration  isolation  of  the  SSA,  and  the  room  temperature  shaker  used 
to  perform  the  calibration  of  the  SSA. 

3.3.1  Thermal  and  Mechanical  Design  of  the  Cryostat 

The  SSA  is  operated  in  a  very  quiet  cryogenic  vacuum  space.  A  low  boil-oflf 
super-insulated  liquid  He  dewar  (William  G.  Goree,  Inc.,  Los  Altos,  California)  was 
constructed  with  a  careful  design.  The  dewar  has  an  inner  diameter  of  61  cm  and  an 
inner  depth  of  1.8  m.  This  large  cryogenic  space  is  carefully  shielded.  The  inner  jacket 
of  the  dewar  is  constructed  from  thin  aluminum  and  Fiberglas.  The  upper  section  of 
the  jacket  is  made  of  thin  Fiberglas  which  has  a  low  thermal  conductivity  to  mechanical 
strength  ratio.  The  lower  section  is  made  of  aluminum.  This  combination  helps  to  stiffen 
and  electrically  shield  the  lower  section  of  the  dewar  without  sacrificing  the  hold- time 
of  the  dewar.  A  long  hold-time  of  liquid  helium  is  needed  for  uninterrupted  operation 
of  the  SSA  and  the  SGG.  This  is  also  important  in  reducing  the  boiling  noise  from  the 
liquid  helium. 

A  super-insulated  dewar  was  chosen  over  a  more  traditional  nitrogen  jacketed 
dewar  in  order  to  remove  the  noise  of  the  boiling  liquid  nitrogen.  This  noise  is  of 
much  more  significance  than  the  boiling  helium  since  the  liquid  nitrogen  would  form  the 
primary  heat  sink  for  the  thermal  radiation  in  the  laboratory.  With  the  current  setup, 
no  correlation  between  the  activity  in  the  helium  bath  and  the  noise  in  the  SSA  has  been 
seen. 

The  outer  jacket  of  the  dewar  is  constructed  from  aluminum,  and  the  entire  unit 
is  extremely  rigid.  This  rigidity  is  important  in  order  to  reduce  the  long  term  “angular" 
drift  of  the  dewar.  Electromagnetic  disturbances  are  also  an  important  source  of  noise 
in  a  sensitive  experiment.  The  cryostat  has  been  carefully  designed  with  this  in  mind. 
In  addition  to  outer  and  inner  layers  of  aluminum,  three  addiliona!  layers  of  shielding  lie 
between  the  inner  and  outer  jacket  of  the  dewar.  Two  layers  of  mu-metal  insure  that  any 
external  magnetic  field  from  the  earth  will  bo  almost  entirely  excluded  from  the  cryogenic 
space.  A  third  shield,  constructed  from  a  thin  layer  of  lead  foil  in  close  contact  wdth  the 
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inner  aluminum  jacket  of  the  dewar,  forms  a  superconducting  shield  which  will  exclude 
any  residual  electric  and  magnetic  field.  The  resulting  cryogenic  space  is  exceptionally 
free  of  electromagnetic  disturbances. 

The  cryostat  insert  which  fits  snugly  into  the  cryogenic  space  is  shown  in  Fig. 
3.9.  The  top  plate  of  the  insert  bolts  to  the  top  of  the  cryostat  and  a  large  rubber 
0-ring  seals  the  cr_yogcnic  space  from  the  laboratory  atmosphere.  This  allows  the  liquid 
nitrogen,  which  is  initially  transferred  into  the  dewar  to  precool  the  insert,  to  be  pumped 
on  by  an  external  roughing  pump  and  brought  close  to  its  triple  point  (63.15°K).  This 
lowers  the  temperature  of  the  insert  and  makes  for  a  more  efficient  transfer  of  helium. 

The  large  aluminum  vacuum  can  at  the  base  of  the  insert  is  anodized  and  coated 
internally  with  a  thin  layer  of  lead.  This  lead  is  evaporated  onto  the  inner  surface  of  the 
vacuum  can  by  passing  a  large  current  though  a  niobium  boat  which  is  filled  with  lead. 
This  lead  thin  film  forms  a  type-I  superconducting  shield  which  further  reduces  magnetic 
flux  level  inside  the  experimental  space.  A  final  coating  of  Kodak  KPR  photo-resist  is 
added  to  protect  the  thin  coating. 

Seven  hollow  Fiberglas  tubes  support  the  vacuum  can  from  the  brass  plate  at 
the  top  of  the  insert.  Six  of  these  tubes  are  equally  spaced  in  a  ring  near  the  edge  of  the 
vacuum  can  and  the  seventh  lies  at  the  center.  Each  end  of  the  hollow  Fiberglas  tubes  is 
fitted  to  an  aluminum  fixture  which  is  used  to  bolt  the  tubes  in  place.  These  aluminum 
fixtures  are  vacuum  sealed  with  a  thin  layer  of  Armstrong  A 12  epoxy  (Lunax  Products, 
Inc.,  Fullerton,  California),  which  matches  the  thermal  expansion  between  the  Fiberglas 
and  the  aluminum.  These  Fiberglas  tubes  in  turn  support  several  radiation  shields 
constructed  from  copper  coated  Fiberglas  PC- board  material.  These  radiation  shields 
provide  thermal  shielding  and  mechanically  stiffen  the  insert.  The  spacing  between  the 
radiation  shields  is  largest  just  above  the  vacuum  can  and  decreases  towards  the  top  of 
the  cryostat,  where  the  thermal  gradient  is  greatest. 

All  the  leads  that  enter  the  cryostat  are  vapor  cooled  and  heat  sunk  to  the 
liquid  helium  bath.  Two  types  of  leads  enter  the  cryostat:  high  current  leads  which 
are  constructed  from  copper  magnet  wire,  and  low  current  leads  which  are  constructed 
from  manganin  wire.  Copper  has  a  much  higher  thermal  conductivity  than  manganin. 
For  this  reason,  copper  was  used  only  for  those  few  leads  that  must  carry  high  current. 
All  leads  were  spirally  wound  around  the  Fiberglas  support  tubes.  This  increases  the 
surface  area  of  the  wire  and  allows  for  greater  vapor  cooling  to  take  place.  In  addition, 
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the  copper  leads  were  heat  sunk  to  the  radiation  shields. 

The  dewar  is  equipped  with  a  helium  level  gauge  (American  Magnetics,  Inc., 
Oak  Ridge,  Tennessee),  which  monitors  the  level  of  helium  above  the  top  of  the  vacuum 
can.  In  addition,  several  carbon  resistors  with  large  temperature  coefficients  are  used  to 
monitor  the  level  of  helium  or  nitrogen  in  the  cryogenic  space.  A  large  nitrogen  boil-off 
heater  suspended  below  the  vacuum  can  is  used  to  boil  away  the  liquid  nitrogen  prior  to 
the  initial  filling  with  liquid  helium.  The  heater  is  a  thin  aluminum  plate  with  thirty  5 
W,  1  fi  resistors  epoxied  to  its  surface  and  connected  in  series. 


3.3.2  Electrical  Isolation  of  the  SSA 

A  large  twelve-sided  aluminum  box  on  top  of  the  dewar  insert  contains  filters 
for  all  leads  entering  the  cryogenic  space.  Each  face  of  the  filter  box  (except  one)  has 
one  hermetically  sealed  forty-one  pin  connector  bolted  into  it.  The  one  face  without  a 
connector  contains  the  vent  line  for  the  dewar.  A  twelve-sided  copper  wall  fits  concen- 
tricly  inside  the  box.  Each  face  of  the  copper  wall  (except  the  one  with  the  vent  line) 
has  forty-one  coaxial,  bushing  mounted  emi  filters  (Model  9001-100-1010,  Murata-Erie 
Corp.,  Toronto,  Ontario)  passing  through  it.  These  filters  are  high  quality  pi  filters  con¬ 
taining  two  capacitors  and  one  inductor.  Above  thirty  to  forty  kllz,  each  filter  reduces 
the  rf  noise  by  80  dB.  Presently,  only  five  of  the  eleven  connectors  are  being  used.  The 
remaining  six  will  be  used  by  the  SGG  when  it  is  integrated  with  the  SSA.  The  cen¬ 
ter  of  the  filter  box  is  hollow  and  houses  room  temperature  accelerometers  used  in  the 
calibration  of  the  SSA. 

In  addition  to  the  rf  filters  at  room  temperature,  there  are  eleven  PC  boards 
which  contain  filter  capacitors.  These  are  mounted  on  top  of  the  vacuum  can  in  the 
helium  space.  These  capacitors  provide  additional  filtering  and  help  insure  that  the  rf 
level  of  noise  in  the  vacuum  sp^lce  is  extremely  low. 

Each  pair  of  feedback  and  oscillator  leads  that  connect  to  the  SSA  are  fed 
through  a  separate  hollow  Monel  tube  that  extends  from  the  brass  plate  at  the  top  of 
the  dewar  insert  to  the  top  of  the  vacuum  can.  In  this  way  the  twin-ax  leads  are  shielded 
all  the  way  to  the  top  of  the  vacuum  can. 
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3.3.3  Vibration  Isolation  of  the  SSA  and  Shaker  Design 


A  room  temperature  shaker,  which  is  free  to  translate  in  x,  y  and  z,  and  rotate 
about  the  z  axis,  is  used  in  the  calibration  and  test  of  the  SSA.  A  cross  section  of  this 
shaker  is  shown  in  Fig.  3.10.  The  shaker  is  mounted  on  top  of  the  rf  filter  box  on  the 
cryostat  insert.  The  center  mass  of  the  shaker  is  suspended  in  vacuum  by  a  spring  from 
the  top  of  a  Plexiglas  tower.  The  center  mass  of  the  shaker  is  machined  in  the  shape  of  a 
cross.  A  strong  rare  earth  magnet  is  mounted  at  the  end  of  each  of  the  four  protruding 
arms  of  the  cross.  Each  of  these  magnets  sits  between  two  coils  which  are  mounted  to 
the  outside  of  the  shaker  housing.  By  driving  current  through  these  eight  coils  in  the 
appropriate  combinations,  the  central  mass  of  the  shaker  can  be  driven  in  three  degrees 
of  freedom:  linear  in  x  and  y,  and  angular  about  the  z  axis  (vertical).  One  more  magnet 
mounted  above  the  cross  and  below  the  supporting  spring  is  used,  along  with  a  set  of 
matching  coils  mounted  outside  the  Plexiglas  tower,  to  drive  the  shaker  along  the  ;  a.xis. 
The  housing  of  the  shaker  is  constructed  entirely  of  Plexiglas.  This  simplifies  the  task 
of  alignment,  and  allows  us  to  measure  the  rotational  displacement  and  calculate  the 
angular  acceleration  applied  by  the  shaker  to  tb<^  SSA. 

The  calibration  of  acceleration  applied  to  the  SSA  by  the  shaker  in  the  lin¬ 
ear  degrees  of  freedom  is  provided  by  a  set  of  three  room  temperature  accelerometers 
(Sunstrand  Data  Corp.,  Redmond,  \VA)  mounted  immediately  below  the  shaker  mass. 
A  thin  Fiberglas  rod  is  attached  to  the  set  of  accelerometers  and  travels  down  through 
the  hollow  Fiberglas  rod  at  the  center  of  the  dewar  insert.  This  Fiberglas  rod  extends 
approximately  12  cm  below  the  top  of  the  vacuum  can  and  is  used  to  support  the  SSA. 

The  combined  structure  which  supports  the  SS.\  can  l)e  modeled  as  a  compound 
pendulum  which  is  supported  by  a  spring,  i'he  mode  slrnrture  of  this  support  is  in 
general  quite  complicated,  but  in  practice,  we  are  only  concerned  with  the  frequency 
response  well  below  the  resonance  modes  of  the  suspension  structure.  At  these  low 
frequencies,  the  modes  of  the  suspension  structure  can  be  ignored.  .4t  higher  frequencies, 
the  suspension  can  be  thought  of  as  a  second  order  low-pass  fdter.  This  filtering  of  high 
frequency  noise  is  important  in  order  that  the  dynamic  range  of  llie  S.S.\  not  be  exceeded 
and  the  SQUID  unlocked. 


P’ignrf'  .1.10.  ;\  cross  sectional  view  of  tlic  S.SA  shaker. 
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3.4  Computer  Control  and  Interface  Electronics 


The  SSA  contains  nine  pairs  of  current  charging  leads,  sixteen  heat  switches, 
six  pairs  of  current  sensing  leads,  and  six  pairs  of  feedback  control  leads.  Given  the 
large  number  and  combinations  of  parameters  to  control  simultaneously,  much  of  tlie 
experimental  control  htid  to  be  automated. 

An  AT-compatible  computer  forms  the  central  core  of  the  control  system  used 
with  the  SSA.  This  computer  is  equipped  with  two  four-channel,  sixteen-bit  data  ac¬ 
quisition  boards  (Model  DT2816,  Data  Translation,  Inc.,  Marlboro,  MA)  which  allow 
the  simultaneous  capture  of  the  six  outputs  from  the  SSA.  A  general  purpose  GPIB 
card  (Model  GPIB-PCIII,  National  Instruments  Corp.,  Austin,  TX)  is  used  to  control  a 
programmable  current  supply  (Model  2500EP,  Valhalla  Instruments,  San  Diego.  C.4),  a 
precision  DVM  (Model  195A.  Keithley  Instruments,  Inc.,  Cleveland,  OH),  and  six  lock-in 
amplifiers  (Model  3961B,  Ithaco  Scientific  Instruments  Inc.,  Ithaca,  NY).  A  parallel  I/O 
card  (Model  DT2817,  Data  Translation,  Inc.,  Marlboro,  MA)  is  used  to  control  custom 
cryostat  interface  electronics. 

A  custom  interface  control  containing  switchable  current  supplies,  high  and 
low  current  relays,  and  two  separate  instrumentation  amplifiers,  was  constructed  by  the 
Electronics  Development  Group  at  the  University  of  Maryland.  This  interface  box  is 
controlled  manually  or  by  twenty-four  optically  isolated  digital  logic  lines  connected  to 
the  AT-compatible  computer. 

The  front  panel  of  the  interface  can  accept  the  input  from  two  separate  pro¬ 
grammable  current  supplies.  Presently,  only  one  of  these  is  in  use,  but  changes  in  the 
design  of  the  SSA  in  the  future  may  necessitate  the  use  of  two  separate  supplies.  High 
current  relays  are  then  used  to  route  current  to  the  appro[)riale  charging  leads  of  the 
SSA.  Low  current  relays  are  then  used  to  control  which  heat  .switches  in  the  SS.\  are 
pulsed  on  and  off.  The  amount  of  current  used  to  drive  each  of  these  heat  switches  is 
set  manually  at  the  front  panel.  Two  separate  instrumentation  amplifiers  (one  for  each 
of  the  separate  programmable  current  supplies)  are  u.sed  to  differentially  amplify  the 
voltage  across  the  current  charging  leads  while  current  is  being  stored  in  the  SSA.  These 
outputs  are  used  to  control  current  protection  circuits  which  shunt  current  away  from 
the  SSA  if  a  preset  voltage  level  is  exceed. 


71 


Chapter  4 

Test  of  the  Superconducting 
Six- Axis  Accelerometer 


In  the  previous  two  chapters  we  described  the  theory,  design,  and  construction 
of  the  SSA.  In  this  chapter  we  will  describe  the  lest  auu  performance  of  the  SSA. 

Initially  the  proof  mass  rests  on  the  three  lower  coil  forms  and  all  six  inductance 
bridges  are  maximally  unbalanced.  Our  task  is  to  levitate  the  proof  mass  and  balance 
all  six  bridges  simultaneously.  If  each  inductance  in  the  SSA  were  identical,  this  task 
would  be  greatly  simplified:  all  levitation  currents  would  be  identical  and  there  would 
be  no  need  to  adjust  the  angular  orientation  of  the  proof  mass.  Every  attempt  has  been 
made  to  insure  that  all  of  the  coils  in  the  SSA  are  as  similar  as  possible.  Despite  this, 
variations  in  the  diameter  of  the  superconducting  wire,  errors  in  the  orthogonality  of  the 
windings,  and  mechanical  alignment  errors  conspire  to  insure  that  each  coil  is  unique. 
This  uniqueness  forces  all  three  levitation  currents  to  be  slightly  different  and  additional 
currents  must  be  stored  in  each  of  the  angular  levitation  circuits  in  order  to  balance  the 
three  angular  sensing  circuits 

This  balancing  process  is  straightforward,  but  it  must  he  carried  out  in  an 
iterative  manner.  Each  time  one  of  the  currents  in  the  SS.A  is  adjusted  in  order  to 
balance  one  of  the  sensing  bridges,  the  position  of  the  proof  mass  changes.  This  small 
change  in  position  causes  the  current  stored  in  each  of  the  other  five  levitation  circuits 
to  change.  This  change  shifts  the  position  in  tho.se  other  five  degrees  of  freedom,  and 
forces  us  to  store  new  currents  in  each  of  those  five  levitation  circuits. 
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In  addition  to  the  control  for  the  SSA,  the  interface  box  contains  additional 
circuitry  to  monitor  the  conditions  inside  the  cryostat.  Additional  high  current  relays 
control  the  nitrogen  boil-off  heater  and  allow  the  computer  to  control  the  nitrogen  boil-off 
cycle. 

The  SSA  instrument  has  been  designed  to  be  completely  automated.  Presently, 
all  currents  in  the  SSA  are  stored  under  computer  control,  current  protection  circuitry 
can  be  checked  and  reset  under  computer  control,  and  boil-off  of  the  nitrogen  in  the 
dewar  is  automated. 

Each  lock-in  amplifier  contains  a  precision  oscillator  which  is  amplitude  and 
frequency  controllable.  These  oscillators  are  used  to  drive  the  sensing  circuit  in  the  SSA. 
Presently,  the  frequency  and  amplitude  level  of  each  oscillator  is  set  manually  and  the 
complete  automation  of  this  is  underway. 
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This  process  could  go  on  indefinitely,  were  it  not  for  the  orthogonal  design  of 
the  SSA.  Each  time  one  of  the  sensing  bridges  in  rebalanced,  the  change  in  position 
along  that  degree  of  freedom  is  much  greater  than  the  change  in  position  in  the  other 
five  orthogonal  directions.  In  practice,  iteration  through  four  or  five  cycles  often  brings 
us  to  the  point  where  we  are  limited  by  the  resolution  of  the  current  supply  used  to  store 
currents  in  the  SSA. 


4.1  Accelerometer  Parameters 

The  experimental  parameters  for  the  SSA  are  summarized  in  Table  4.1.  The 
masses  of  the  proof  mass  and  SSA  assembly  are  measured.  The  moments  of  inertia  are 
computed  from  the  mechanical  design.  The  levitation  and  sensing  inductance  parameters 
are  calculated  in  Section  4.1.1.  The  inductance  of  the  feedback  transformer  secondary, 
L,  was  measured,  as  described  in  Section  4.1.2,  and  is  an  average  value  typical  of  the 
six  feedback  transformers.  The  input  inductance  of  the  SQUID,  Lsq,  is  taken  from  the 
manufacturer’s  data  sheet.  The  remaining  three  parameters,  c,  di,  and  ds,  are  taken 
from  the  mechanical  design. 

4.1.1  Levitation  and  Sensing  Inductance  Modeling 

From  our  knowledge  of  the  geometry  of  the  superconducting  coils  in  the  SSA, 
we  can  calculate  the  parameters  which  describe  the  inductance  of  each  of  the  coils.  Each 
of  the  circular  windings  in  the  sensing  or  levitation  coil  can  be  thought  of  as  a  closed 
circular  loop  containing  a  current  /.  This  loop  creates  a  corresponding  image  current 
loop  in  the  superconducting  proof  mass.  By  calculating  the  force  between  all  current 
loops  and  image  current  loops  and  summing,  we  can  calculate  the  force  on  the  proof 
mass  from  each  coil.  This  force  will  be  a  function  of  the  distance,  d,  and  orientation,  6, 
of  the  coil  with  respect  to  the  proof  mass.  The  force  acting  between  the  proof  mass  and 
one  of  the  sensing  or  levitation  coils  can  be  related  to  the  variation  in  inductance  with 
position.  By  compari.ig  these  two  forces,  we  can  deduce  the  variation  of  inductance  with 
position. 

The  force  between  *wo  loops  carrying  equal  and  opposite  currents  /,  is  [8] 

r  fio  ,2  I  I  ■  dVi  ,,,, 
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Parameter 

Symbol 

Value 

Mass  of  the  Pcoof  Mass: 

m 

0.144  kg 

Mass  of  the  SSA  Assembly; 

Af 

■L8  kg 

Moir  ent  of  Inertia,  Proof  Mass: 

/ 

4.0  X  10“^  kg  nr 

Moment  of  Inertia,  SSA  Assembly; 

hi 

8.2  X  10-"'  kg  m‘ 

Secondary  Inductance  of  Feedback  Transformers: 

L 

240  fiU 

Inductance  of  SQUID  input  Coil: 

2.0  ^i\l 

Parallel  Displacement  of  Inductance  Coil  Center 
from  Proof  Mass  Center; 

c 

1-3.5  X  10~^  m 

Perpendicular  Displacement  of  Inductance  Coil  Center 
from  Proof  Mass  Center 

Levitation  Coil: 

dp 

1.69  X  10-'*  in 

Sensing  Coil: 

d,s' 

1..51  X  10-'*  rn 

Levitation  Coil  Inductance  Parameters; 

Ll. 

31  /til 

A/, 

6.1  X  iO"^  ll/m 

■h. 

1.8  X  10-^  H 

48  H/rn'^ 

Sensing  Coil  Indue^ance  Parameters: 

Ls 

1.8  /rll 

A.s 

5.2  X  10-^  11/in 

■h 

2.1  X  lO"'*  H 

1  S 

5.9  11  /nr 

T;.'  Ir  l.l.  IOxixt' metal  parameters  for  tlie  SSA. 


where  rfl,  and  (i\<  ari’  infiniteiitiiai  S('Km"ut^^  of  loop,->  1  and  re.vpei  tivf  ; y,  and  Xi.  i' 
the  vector  joining,  tho>e  two  M'enicuts.  Our  two  current  loops,  one  pliysica!  of  radiius  /f j 
and  one  image  of  radius  /O-  are  >hown  in  fig.  1.1.  hach  loop  situated  at  a  di.sia:iie 
d  +  s  and  orientation  d  with  respect  to  the  supercotiduct  ing  jdane.  .Xlthoiigh  If.j.  i  1 . !  ' 
contains  force  components  in  the  ly  and  directioii.s.  tlie'-<>  (om[)o,nents  must  ^.iiiish, 
I  herefore,  we  need  to  <  oinpute  only  tlie  x  cotnponenl  oi  tie'  I'ni.-.  lApanding  1  (].  '  I!  ‘ 
t  h  rough  h  r-t  order  in  x  and  d.  and  sii  mini  tig  over  all  loop-  in  hot  h  -u  jier  vi''id  u,  t  i  ng  roi  I  - . 
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Figure  4.1.  A  schematic  picture  of  a  current  loop  above  a  superconducting 
plane  interacting  with  an  image  current  loop  below. 


we  find 


^  /  _  i?i  /?2  COs(^2  ~  ^1  ) _ 

I  [.i(P  +  R2+  j^2  _2R,R2  C0S(^2  -  01  )] 


X  \  2d  2x  I  I  — 


[A(P-  ^  R]^  R\-  2R^R2Cos(92  -0v)]\ 


d9id92,  (4.2) 


where  we  have  dropped  all  terms  which  vanish  due  to  symmetries  in  the  integrals. 
The  torcpie  apjilied  to  loop  1  In-  loop  2  is  given  by  the  expression 


,  /'o  , 

12  =  — / 
4;r 


'// 


R]  X  Xi2d]i  ■  dlj 


where  Ri  is  the  vector  from  the  center  of  loop  1  to  the  infinitesimal  line  segment  r/l]. 
An  expansion  of  this  exi)ression  for  small  x  and  9.  and  summation  over  all  loops  .  'ach 
coil,  gives 


n'l  //, 


ros(^2  ■■  9\)  ros9id9idt)2 
[  irP  /^f  +  R]  -  2/f,  R2  ros(t>2  -  ^1 


R  ^  cos  9i  +  /(‘jcosf'^) 

[  |7t  r  /f(’ 7  Rj  Ylfjt.  cos(fl.“^  f?,  )j 


-  /7-  cos  9,  >  . 


The  forces  acting  on  the  proof  mass  can  also  be  calculated  from  the  variation  of  induc¬ 
tance,  //,  with  position: 


fr 

n 


1  2^ 
2  dx' 
\  .dL 

2  de' 


(4.5) 

(4.6) 


Substituting  the  expression,  Eq.(2.89),  for  the  inductance  into  these  two  equation.s,  we 
find 


F.  -  ^A/2-i7x/2,  (4.7) 


Setting  Eqs.  (4.5)  and  (4.6)  equal  to  Eqs.  (4.2)  and  (4.4),  and  solving  for  and  7, 
we  find 

R\  /?2  cos(02  - 


A 


= 


TT  J  J 


7  = 


'‘“I'T'// 

R 

u 


TT  J  J 

Ri  R2 


[4(P  +  Rl  +  Rl-  27?,  7?2  cos(02  - 

7?,  7?2COs(02  -  ) 

[4(P  +  7?2  +  7?|  -  27?,  7?2  cos(6'2  - 


(iB\  (i&2 


(4.9) 


I  [4(P  +  7?2  +  7?2  -  27?i7?2Cos(02  -  )] 

^  R'^R'y  COS(»2  -  ^1)  cos  0,(7^,  (/02 

^  n  n  J  J 


(4.10) 


fl.  R,  —  [4(72+;?2+fi2_.27?,7?2COS(02-^l)]"^' 
6(7?,  cos  0,  4- 7?2  cos  02)d^  ^ 


\  [4d2  +  7?2  +  7?2  -  27?,7?2Cos(^2  -  Pi)]  j 

The  expressions  for  A  and  7  can  be  reduced  to  one-dimensional  integrals 

7?,  7?2u  cos  u 


ft,  ftj 


7  =  '^Ro 


[4(72  +  7^2  +  7?2  -  2/?,  /?2  cos 
7?)  7?2  cos  u 


--  d  n , 


(4.11) 


(4.12) 


ft,  ft2  [4d^  +  /??  +  /?.2  -  2/?,  7?2  cos  u 
f  !2(7'2 


3/2 


[  [4(7'2  +  7?2  +  /?2  _  27?,  7?2  cos  !/] 


-  1  >  du. 


(4.13) 


and  numerically  integrated  using  Romberg's  me' hod  of  order  10  [13]. 

The  expression  for  /7  is  less  tractable  and  must  be  handlc'd  a.s  a  two-dimensional 
integral.  This  was  implemented  mumerically  as  a  one-dimensional  integrid  which  per 
forms  a  series  of  one-dimensional  integrals.  'I'hese  routines  wore  implemented  in  the  (' 
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language,  and  the  calculation  of  /?  for  the  levitation  coil  took  in  excess  of  24  hours  of 
CPU  time  on  a  VAX  8800. 

Each  sensing  coil  consists  of  30  turns  of  Nb-Ti  wire.  The  inner  and  outer 
diameter  of  each  sensing  coil  is  0.51  cm  and  1.22  cm,  respectively.  Once  the  proof  mass 
is  levitated,  the  center  of  each  wire  in  the  sensing  coil  rests  330  /im  away  fiom  the  proof 
mass  surface.  Substituting  these  parameters  into  our  expressions  for  A,  0,  and  7,  we 
find,  after  numerically  integrating. 


As 

=  5.4  X  10"^  ll/m. 

(4.14) 

0s 

=  2.1  X  10“'’  H, 

(4.15) 

~)S 

=  5.9  ll/m^. 

(4.16) 

.Multiplying  A5  by  the  sensing  coil  to  proof  mass  separation,  330 //m,  gives  Ls  =  1.8 /rll. 

Each  levitation  coil  consists  of  two  layers  of  Nb-Ti  wire.  Each  layer  holds  40 
turns.  The  inner  and  outer  diameter  of  each  levitation  coil  are  1.28  cm  and  2.11  cm. 
respectively.  Once  the  proof  mass  is  levitated,  the  centers  of  the  wires  which  make  up 
the  upper  layer  in  each  levitation  coil  rest  at  a  distance  dj  -  440  /zm  away  from  the 
proof  mass  surface,  and  the  centers  of  the  wdres  which  make  up  the  lower  layer  rest  at 
a  distance  ^2  =  570  /zm  away  from  the  proof  mass  surface.  Because  each  levitation  coil 
consists  of  two  layers,  the  expressions  for  A,  7,  and  0  are  slightly  more  complicated  and 
the  term  4d^,  which  ropre.sents  the  distance  between  the  physical  current  and  the  imago 
current,  must  be  replaced  by  4f/i,  4^21  and  (dj  +  ^2)^  as  the  double  summation  runs  over 
the  various  layers  of  physical  current  and  image  current.  Expanding  the  expressions  for 
A,  7,  and  0  into  these  three  cases,  and  substituting  the  geometrical  parameters,  we  find, 
after  numerically  integrating. 


Al 

=  6.1  X  10" 

^  H/m, 

(4.17 

0L 

=  1.8  X  10" 

Ml. 

(4.18 

11. 

=  48  1I/in^ 

(4.19 

Multiplying  A/^  by  the  average  levitation  coil  proof  nias>.  separation,  (d|  +  give.s 

il  =  31  /zll. 

Each  of  the  expressions  for  .\,  7,  and  3  ignores  the  rearrangz'inent  of  the  current 
in  each  of  the  siipturc,  .dnrt ing  wires  due  the  oIIhu  su [X'rronduct ing  wires  and  image 
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currents,  and  neglects  the  finite  size  of  the  superconducting  plane.  These  are  small 
errors,  and  our  estimation  of  the  coil  parameters  (  Psi  A5,  75,  Pi,  A/,,  and  7^)  should 
be  correct  within  ten  to  twenty  percent. 

4.1.2  DC  Inductance  Measurements 

By  storing  a  small  current  in  each  of  the  levitation  circuits  and  then  allowing 
that  current  which  is  trapped  in  the  circuit  to  decay,  it  is  possible  to  measure  the  amount 
of  flux  trapped  in  each  of  the  superconducting  loops.  Since  the  amount  of  flux  trapped 
in  a  superconducting  loop  is  the  product  of  the  inductance  and  the  current  stored,  it  is 
possible  to  deduce  the  circuit  inductance. 

This  technique  [lUJ  integrates  the  voltage  generated  across  the  heat  switch,  i.e.. 

f+00  y  +  oo 

$  =  Ltotf  =  J  V{t)dt  ~  -  j  (4.20) 

where  Ltot  is  the  total  inductance  of  the  superconducting  loop.  Ltot  for  each  of  the  three 
linear  inductance  circuits  can  be  computed  from  Eqs.  (2.92)  thru  (2.115)  and  Fig.  2.5. 
We  find 

I'tOtrj:  ~  ^'totry  =  ^'totr^  =  (d-^l) 

Storing  small  currents  in  each  of  the  three  linear  levitation  circuits  docs  not  levitate  the 
proof  mass,  and  Li  in  Eq.  (4.21)  must  be  computed  using  the  spacing  (250  /im)  between 
the  proof  mass  and  the  levitation  coil  when  the  proof  mass  is  resting  on  the  sensing  coil. 
We  find  Li  =  Ay,  250 /tm  =  16//H.  Experimentally,  Ltot,^  =  330  /iH,  =  280  yrll, 

and  Ltotr,  =  300  /iH.  Substituting  the.se  values  into  the  above  equation,  solving  for  L 
and  averaging,  we  find 

/.  =  240  /1 H.  (4.22) 

Flux  measurements  of  the  three  angular  levitation  circuits  are  consistent  with  this  n'sult. 

4.2  Accelerometer  Resonant  Modes 

The  final  set  of  currents  stored  in  the  SS.A,  which  balance  all  six  inductance 
bridges  simultaneously,  are  given  in  Table  4.2.  All  three  linear  levitation  currents  are 
slightly  different  from  each  ^ther,  as  expected.  Initially,  we  attempted  to  store  0.9  \  in 
all  three  angular  levitation  circuits  and  we  discovered  that  one  of  the  heat  switches  in 


the  6^  levll.itiot!  rirr^i'  nc*  ■p-'-ra’ioi  -U.  1  !ii»  jorccd  'if-  -o  •Hye  d  s^'Tririi- i  :it:  pair  of 
levuauo..  .....  '  -  I- 

From  Eqs. (2.118),  (2.119),  (2.120),  (2.124),  (2.12-3),  ana  U.izo;,  .1.^ 
levitation  force  is 


fOCr^  =  1 

Al, 

(4.2.3) 

foCr,  = 

K  -  uk  +  ^k)_ 

Al, 

(4.24) 

fOCr.  =  1 

-  (Ik  +  fk\ 

Al. 

(4.25) 

Setting  these  equations  equal  to  the  applied  external  force,  mgEf^/Z,  substituting  in  our 
experimental  currents,  and  solving  for  A/,,  we  find 

Ac  =  6.6xlO“^H/m  (from/cc,^),  (4.26) 

Al  =  6.9xlO“^H/m  (from/£)Crj,)>  (4-27) 

Al  =  6.8xlO~^H/m  (from/oCrJ-  (4.28) 

The  first  of  these  experimental  values  is  slightly  lower  than  the  other  two.  This  can  be 
understood  if  we  consider  the  fact  that  the  diameter  of  the  Nb-Ti  wire  varies  slightly 
between  coils,  which  changes  the  density  of  the  windings  and  modifies  A^,.  The  average 
experimental  value,  Ac  =  6.8  x  10“^  H/m,  is  11%  higher  than  our  theoretical  value,  but 
is  reasonable  when  we  take  into  account  the  crudeness  of  our  inductance  model. 

The  spring  constants  for  each  of  the  six  degrees  of  freedom  can  be  obtained 
from  the  levitation  potentials  in  Eqs.  (2.118),  (2.119),  (2.120),  (2.124),  (2.12-5),  and 
(2.126).  We  find 

fcr,  = 


kr.  = 

= 


SAl 


4Le  +  7/ 


+  7L  +  2(/26I  +  ifg  ) 


2  y  ,  IT  ,  /  \2 

i  i  29,  +  19,)  ^ 


SAl 


^2 


(4.29) 

A  2  j  *2 

(4.30) 

(4.31) 

[A[^d[,  +  {2ll^  +  2/^,  +  +  f^o,  +  ) 
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Parameter 

Symbol 

Value 

X  Levitation  Current 

hx 

2.639  A 

y  Levitation  Current 

2.. 588  A 

z  Levitation  Current 

hx 

2.455  A 

6x  Levitation  Circuit  1  Current 

hOx 

-0.280  A 

6x  Levitation  Circuit  2  Current 

hox 

0.280  A 

&y  Levitation  Circuit  1  Current 

hOy 

0.786  A 

Oy  Levitation  Circuit  2  Current 

hOy 

1.014  A 

Ox  Levitation  Circuit  1  Current 

he. 

0.918  A 

Ox  Levitation  Circuit  2  Current 

he. 

0.882  A 

X  Resonance  Frequency 

fr. 

21.6  Hz 

y  Resonance  Frequency 

21.5  Hz 

z  Resonance  Frequency 

fr. 

18.2  Hz 

Ox  Resonance  Frequency 

fe. 

14.96  Hz 

Oy  Resonance  Frequency 

fOy 

15.15  Hz 

Ox  Resonance  Frequency 

fe. 

16.23  Hz 

Mechanical  Quality  Factor  for  the  Linear  Degrees  of  Freedom 

Qr 

220 

Mechanical  Quality  Factor  for  the  Angular  Degrees  of  Freedom 

Qe 

700 

Table  4.2.  Levitation  currents,  re.sonance  frequencies,  and  mechanical 
quality  factors  for  each  of  the  six  inodes  of  the  SSA. 


cWl 


L\jc 


2 


(2/^,  +  2/^_  +  l2e,  +  ^ 

^2*2  r  A  2  2 

+  2{lle  +  ^ 


10V 


+  (hSy  -  fl0„) 


{J^iAl  +  +  /^/,)  (2^rr  +  +  ^?9j,  +  ^20,  +  ) 


+  2(/2e^  +  ) 


2  A  2 

ii.  +  i 


+  (ho,  -  ho,) 


I-\IA 

^■l(1^i.  +  L ) 


(4.32) 


(1..33) 


(4.34) 


Substituting  in  our  experimental  currents,  the  average  experimental  value  of  A; .  and 
our  theoretical  parameters,  and  computing  the  resuliant  resonanci'  frequencies,  we  find 


=  21.2  11/.,  (4.35) 

fr^  20.9  11/,  (4.36) 


8.3 


fr. 

=  17.8  Hz, 

(4.37) 

/fix 

=  13.8  Hz, 

(4.38) 

K 

=  13.9  Hz, 

(4.39) 

fe. 

=  14.5  Hz. 

(4.40) 

The  theoretical  linear  resonance  frequencies  match  the  experimental  resonance 
frequencies  quite  closely.  The  theoretical  angular  resonance  frequencies,  however,  are  all 
approximately  10%  below  the  corresponding  experimental  values  This  implies  that  wo 
have  underestimated  the  angular  spring  constant  by  20%.  This  is  significant,  but  the 
angular  spring  constants  are  much  more  dependent  on  the  nonlinear  inductance  model 
parameters,  as  can  be  seen  from  Eqs.  (4.32)  through  (4.34).  These  parameters  should 
be  more  dependent  on  subtle  effects  that  we  have  neglected  in  our  inductance  modeling. 
The  relative  scaling  between  each  group  matches  quite  well  with  the  experimental  values, 
and  confirms  that  our  basic  understanding  of  the  dynamics  of  the  SS.4  is  correct.  It  is 
also  worth  pointing  out  that  without  the  nonlinear  inductance  coefficients,  f3  and  7,  we 
would  be  unable  to  understand  the  mode  structure  of  the  SSA. 

Thus  far  we  have  neglected  the  contributions  that  the  sensing  currents  make 
to  the  resonance  frequencies  of  the  SSA.  We  will  see  in  the  next  section  that  this  is  an 
excellent  approximation. 

4.3  Calibration  and  Noise  Measurement 

In  this  section,  we  characterize  the  circuitry  used  to  drive  the  six  sensing  circuits 
in  the  SSA,  examine  the  calibration  of  the  SSA,  calculate  the  energy  coupling  coefficients, 
and  compare  the  expected  sensitivity  of  the  SSA  with  the  experimental  results. 

4.3.1  Sensing  Drive  Circuitry 

Each  of  the  sensing  circuits  in  the  SSA  is  driven  by  a  low  noise  current-to- 
voltage  converter  connected  to  an  rf  isolation  transformer.  This  drive  circuitry  is  shown 
in  Fig.  4.2.  A  pair  of  pi  filters  are  placed  in  line  between  the  current-to- voltage  converter 
and  the  rf  isolation  transformer.  The  combined  lead  and  filler  resistance  in  each  leg 
of  the  circuit  is  Rj  =  2.0  and  the  pair  of  pi  filters  are  modeled  as  two  capacitors. 
Cf  =  0.18.5 /if’,  and  two  inductors,  Lj  =  0.17  mil.  The  outputs  from  the  two  pi  filters 
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Figure  4.2.  The  rf  isolation  transformer  and  circuitry  used  to  drive  each 
of  the  sensing  circuits  in  the  SSA. 


connect  to  a  tank  capacitance  C  placed  in  parallel  with  the  isolation  transformer.  The  rf 
isolation  transformer  is  wound  on  a  hollow  aluminum  fixture  and  consists  of  a  primary 
and  secondary  coil  of  inductance  Z-j  and  T2,  respectively.  The  mutual  inductance  of  the 
isolation  transformer  is  Mu-  A  thin  brass  sheet  separates  the  primary  and  secondary 
windings,  and  the  eddy  current  losses  in  this  isolation  layer  are  represented  by  a  resistance 
R,  which  is  placed  between  the  tank  capacitor  and  the  primary  of  the  transformer. 

The  experimental  parameters  for  each  of  the  sensing  circuits  are  shown  in  Table 
4.3.  The  capacitances  listed  are  the  nominal  values,  as  measured  at  room  temperature. 

It  is  straightforward  to  show  that  the  current  in  the  primary  of  the  rf  isolation 
transformer  induces  a  current  I2  in  the  secondary  of  the  Iransfori*...:,  of  magnitude 

The  effective  in[)ut  inductance  of  the  isolation  tr<insformer  can  also  be  com¬ 
puted: 

1.2  +  I.S 

This  effective  inductance  will  resonate  with  the  combined  capacitance,  C  +  r'y,  at  a 
frequency 

.  1  1 


L„ff  = 


-'\Aw/ic  +  o)’ 


Circuit 

Li 

L2 

Ml  2 

C 

l^.nl 

feip 

Tl 

3.75  mil 

0.17  mH 

0.53  mH 

lOfiF 

86 

1.00  kHz 

3.54  mH 

0.19  mH 

0.51  mH 

2  fiF 

95  n 

2.12  kHz 

Tz 

3.54  mH 

0.20  mH 

0.50  mH 

1.3  fiF 

10911 

2.47  kHz 

Or 

3.51  mH 

0.17  mH 

0.51  mH 

45 /iF 

52  fl 

0.50  kHz 

Oy 

3.57  mH 

0.19mH 

0.51  mil 

5  fiF 

90 

1.41  kHz 

Oz 

3.54  mH 

0.20  mH 

O..53  mH 

3  /cF 

88  n 

1.84  kHz 

Table  4.3.  Parameters  for  rf  isolation  circuits  in  the  SSA. 

It  can  also  be  shown  that  the  transfer  function  between  the  input  current  into 
the  rf  isolation  circuit,  /,  and  the  current  passing  through  the  primary  of  the  isolation 
transformer,  /j,  is 


h 


1 


(4.44) 


/  1  —  Lef  f{C  +  Cf)  +  Riu)(C  +  C  j) 

where  we  have  neglected  the  contribution  from  the  first  filter  capacitor,  Cj.  The  first 
filter  capacitor  will  only  play  a  significant  role  in  the  drive  circuit  when  the  frequency 
of  interest  is  greater  than  10  kHz,  and  all  of  our  sensing  circuits  operate  at  frequencies 
below  2.5  kHz.  At  resonance,  this  reduces  to 


/i 


1 


-e// 


R]jC-\-  Cj  ■ 


(4.45) 


Combining  this  with  Eq.  (4.41),  we  find  that  the  magnitude  of  the  transfer 
function  relating  the  drive  current,  /,  to  the  sensing  current,  I2.  is 


1 


L 


eff 


M 


12 


RvC  +  C/L2T 


(4.46) 


It  is  clear  from  this  expression  that  the  transformer  damping  resistance  is  crucial  in 
controlling  the  drive  circuitry  gain.  We  expect  this  resistance  to  increase  as  a  function  of 
frequency,  and  this  will  correspondingly  reduce  the  gain  in  the  rf  isolation  transformer 
as  the  frequeiu  y  of  the  ta.nk  circuit  resonance  increases.  This  is  a  rather  rough  model  of 
a  very  complii  ated  distributed  network.  Given  the  sensitivity  of  the  overall  gain  to  the 
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damping  resistance,  we  should  not  expect  our  theoretical  results  to  agree  precisely  with 
the  experimental  results. 

The  damping  resistance,  R,  can  be  computed  by  measuring  the  magnitude  of 
the  circuit  input  impedence,  Using  straightforward  circuit  analysis,  it  can  be  shown 
that  the  input  impedance  of  the  complete  rf  isolation  circuit  is 


—  2Rj  -f  UluL j  -j- 


R  T  iljLqj j 

+  Q)  +  Riu{C  +  Cj) ' 


(-4.47) 


where  we  have  once  again  neglected  the  contribution  from  the  first  filter  capacitor.  At 
resonance,  this  impedance  reduces  to 


2Rj  + 


Lejj  -  2Lj 


Solving  for  R,  we  find 


(4.48) 


R  = 


(C  +  C/) 


OR 


(4.49) 


Substituting  the  experimental  parameters  from  Table  .1.3  into  our  expressions 
for  R,  Le/j  (Eq.  (4.42)),  resonance  frequency  (Eq.  (4.43)),  and  drive  current  transfer 
function  (Eq.  (4.46)),  wc  find  the  theoretical  parameters  a^  .isted  in  I'able  4.1.  The  the¬ 
oretical  values  for  the  resonance  frequencies  are  all  slightly  arger  than  the  corresponding 
experimental  values.  This  is  not  too  surprising,  considering  that  the  capacitance  value 
is  measured  at  room  temperature,  whereas  the  resonance  is  measured  with  the  circuit  at 
4.2  K. 


4.3.2  Residual  Misbalance 


Each  of  the  sensing  circuits  in  the  SSA  is  balanced  by  adjusting  the  levitation 
current  in  the  six  levitation  circuits  which  surround  the  proof  mass.  .4s  de.scribed  at 
the  beginning  of  this  chapter,  this  process  proceeds  in  an  iterative  manner  until,  after 
several  cycles,  we  are  limited  by  the  resolution  of  the  current  snp[>ly  and  a  small  residual 
misbalance  exists  in  each  of  the  sensing  circuits. 

The  in-phase  component  of  this  misbalance  is  due  to  a  slight  error  in  the  proof 
mass  position,  and  is  inductive  in  nature.  In  practice,  thi.-i  romiJonent  of  the  misbalance 
is  removed  by  sending  a  small  DU  current  into  each  of  the  feedback  transformers  and 
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Circuit 

LefJ 

/theory 

R 

Transfer  Function 

2.11  mH 

1.09  kllz 

2.5  fi 

/r.//=  17.5 

2.18mH 

2.31  kHz 

11  n 

Cz 

2.30  mH 

2.72  kHz 

15  Q 

IrJI  =  6.6 

2.00  mil 

0.53  kHz 

1.04  fi 

2.21  mil 

i.49  kHz 

5fi 

hji=n 

Oz 

2.15  mil 

1.92  kHz 

6.6  n 

hjl  =  8.5 

Table  4.4.  Theoretical  parameters  for  rf  isolation  circuits  in  the  SSA. 

thereby  adjusting  the  position  of  the  proof  maiss  until  each  sensing  bridge  in  the  SSA  is 
nulled. 

The  out-of-phase  component  of  the  misbalance  is  independent  of  proof  mass 
position,  and  is  resistive  in  nature.  To  the  extent  that  each  sensing  coil  in  the  SSA  is 
identical  and  suffers  the  same  resistive  loss,  this  resistive  portion  should  vanish.  In  fact, 
each  of  the  sensing  coils  is  very  well  matched  since  each  coil  is  formed  from  a  single  layer 
of  windings  of  equal  turns.  The  levitation  coils  are  a  two  layer  winding,  the  fabrication 
of  which  is  much  more  difficult  to  control,  and  therefore  these  coils  are  less  well  matched. 

In  the  current  design,  levitation  and  sensing  coils  rest  side  by  side  and  are 
weakly  coupled  together.  Thus,  the  limit  on  how  well  the  sensing  coils  match  comes 
from  the  levitation  coils.  We  are  currently  attempting  to  understand  this  problem,  and 
hope  to  improve  this  in  the  new  design  of  the  SSA. 

In  order  to  quantitatively  understand  the  import  ance  of  the  residual  resistive 
misbaJance,  it  is  interesting  to  consider  the  bridge  circuit  shown  in  Fig.  4. .3.  Three  of  the 
inductances  in  the  bridge  are  identical,  and  suffer  dielectric  and  eddy  current  losses  so 
that  they  present  an  impedence  Rs-  The  fourth  inductor  is  slightly  different  from 

the  other  three,  and  presents  an  impedance  iuLs  +  Rs  +  f^Rs-  Solving  for  the  current 
through  th'.  input  roil  and  the  five  other  bridges,  we  nnd 

hRs  f  cos 

4  [ /?,s'  "t"  i^{()Ls  +  )] 

where  we  hav<  ,'X[)ande(l  through  first  order  in  6Rs- 


(4.50) 
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Measurements  of  the  resistive  component  of  the  impedance  of  a  single  sensing 
coil  at  three  separate  frequencies  give 


/?5(0.5kHz)  =  19..') //P, 

{AM) 

f«.s(1.0kHz)  =  46..3pP, 

(1.52) 

fl5(5.0kHz)  =  171 /iP. 

(4. .53) 

Recalling  that  Ls  -  1.8//H,  we  find,  as  must  be  the  case 

carrier  frequencies  used.  Ap[)lying  this  limit  to  Eq.(‘1.50), 

the  flux  signal  through  the  SQUID  loop  is 

,  lls  <  i^'Ls  for  the  range  of 

we  find  that  the  amplitude  of 

^  Msq  Ills  Ills 

4u)(i)E.s  T  1'Sq)  11  s 

(  l..')l) 

where  Msq  =  2-0  x  10  is  the  mutual  inductance  between  the  SQUID  input  cull  and 
the  SQUID  loop. 

The  slew  rate  of  the  SQUID  feedback  cont  roller  deti'rmines  the  maximum  signal 
size,  4>max,  that  the  SQUID  am()lifier  car.  handh*.  l  lius,  the  maxitnutn  amplitttde  of  the 
current  into  the  .sensing  bridge  is 

^  Iw’IG/y.S'  f  f>.S’y  I'l’ni  iJ  /  1  '  -  'i 

max  - 
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For  the  Quantum  Design  Model  2000  SQUID  controller,  the  maximum  signal 
size  is  8.0  X  10^  $o  up  to  about  60  Hz.  Above  this  it  drops  off  with  frequency,  reaching 
1  at  approximately  20  kHz.  At  2.5  kHz,  the  maximum  signal  is  approximately  20  4»o. 
Assuming  a  5  %  misbalance  in  Rs, 

/mar  =  0.34  A.  (4.56) 

The  actual  limit  on  the  drive  current  through  the  bridge  driven  at  2.5  kHz  was  found  to 
be  close  to  this  value.  Note  that  at  lower  frequencies  ^max  will  be  higher  so  Imax  should 
be  higher.  In  fact,  for  the  bridges  driven  at  the  lower  frequencies,  the  current  was  limited 
by  a  superconducting-to-normal  transition,  probably  caused  by  eddy  current  healing  in 
the  isolation  transformer 

4.3.3  Linear  Acceleration 

The  three  linear  acceleration  circuits  are  calibrated  by  driving  the  shaker  along 
the  vertical  and  comparing  the  signal  from  the  vertically  oriented  room  temperature 
accelerometer  to  the  signals  from  the  SSA.  The  r^,  r^,  and  circuits  were  driven  with 
currents  of  13.25  mA,  0.75  mA,  and  13.25  mA,  respectively.  Multiplying  these  currents 
by  the  three  corresponding  gains  {hjl)  listed  in  Table  4.4,  we  find 


232  mA. 

(4.57) 

5.7  mA, 

(4. .58) 

113  mA. 

(4.59) 

The  Ty  levitation  circuit  was  limited  in  drive  current  because  the  leads  used  to  drive  the  ry 
feedback  transformer  were  damaged  during  the  coo!  down  and  no  feedback  current  could 
be  sent  to  that  circuit.  This  left  a  large  misbalance  signal  which  caused  the  .SQUID  to 
unlock  if  a  current  with  an  amplitude  larger  than  0.75  mA  was  driven  into  the  ry  sensing 
circuit. 


The  expi'riniental 

sensitivities  for 

Cn  ai'd  r-  a1 

:  these  drive  currents  were 

measured  to  he 

3.1  X  lO'’' 

(.1.60) 

9.9  X  10  ’ 

(4.01) 

1.0  :<  lO'’  <I>, ,/,(?. 

(4.02) 
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The  theoretical  transfer  functions  for  the  three  linear  degrees  of  freedom  are, 
from  Eqs.  (2.195)  through  (2.197), 

1 


1 


///  \  =  _ _ _ _ _ _ 


Ir 


A5 


\/2  Lsq  +  6L5  +  iuyi^lQr  -  u.-'^ 

It  A5  1 


(■i.Gd) 

(l.Gl) 

(■1.G5) 


\/2  Lsq  +  GLsi^r  +  ii^r^lQr  - 
These  equations  were  derived  for  identical  resonance  frequencies  and  identical  drive  lev¬ 
els.  All  three  of  our  linear  resonance  frequencies  differ  slightly,  so  that  must  be 
replaced  by  ^rj,>  and  Ur^  in  each  of  the  three  e(iiiations.  All  three  drive  levels  differ, 
and  Ir  must  be  similarly  replaced  by  /rj,,  and  h,-  Making  these  substitutions  and 
taking  the  limit  where  u;  <C  cj^,  we  find 

Irr  As- 


{IIa.,)r 


{Hfict}  rm.T 


4\/2 

Tr  +  ^Ls 

As 

4v/2tr2/ 

7^  ^-s’Q  +  6/>.s' 

Ir, 

A5 

4\/2  TT^y 

■2  Lsq  +  I^Ls 

(4.66) 

(4.67) 

(4.68) 


Substituting  the  values  of  sensing  currents  and  other  parameters  into  these  equations, 
we  find 


=  3.5  X  10^  <ho/(7, 

(4.69) 

=  8.8  X  10^  ^Jg. 

(4.70) 

=  2A  X  10^ 

(4.71) 

The  first  two  of  these  theoretical  sensitivites  agree  witj)  the  experinumtal  values  within 
10%.  This  is  a  rather  good  agreement,  given  the  crudenes.s  with  which  we  are  able  to 
model  the  damping  in  the  rf  isolation  transformer.  The  third  value  is  50%  larger  than 
the  experimental  value. 

A  log-log  plot  showing  the  low-fr<’qnency  noise  spectra  from  0.3  mHz  to  0.3 
Hz  for  the  and  sensing  circuits  is  shown  in  fhg.  4.4.  .Above  0.1  Hz,  the  spectrum 
is  dominated  by  the  background  seismic  noise  of  the  laboratory.  Below  0.1  Hz,  the 
spectrum  appears  to  be  dominated  by  1//  type  noi.se.  The  broad  peak  around  0.2  Hz 
corresponds  to  the  peak  in  the  .seismic  noise  of  the  earth  reported  by  seismologists  [1,2], 
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Figure  4.4.  The  very  low  frequency  spectra  of  noise  as  measured  by  the 
Tj;  and  r,  sensing  circuits. 
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as  is  shown  in  Fig.  4..5.  The  peak  is  attributed  to  the  periodic  pounding  of  the  shores 
by  ocean  waves.  Its  amplitude  is  strongly  dependent  on  a  number  of  factors,  including 
the  distajice  to  the  shore  and  the  amplitude  of  the  ocean  waves,  so  it  is  not  suprising 
that  the  agreement  between  the  peak  heights  on  the  two  figures  is  not  better.  Note  that 
the  seismic  noise  drops  by  over  two  orders  of  magnitude  below  this  peak.  The  fact  that 
the  SSA’s  signal  does  not  reach  a  lower  minimum  is  an  indication  that  the  spectrum  is 
dominated  by  instrument  noise  below  0.1  Hz. 

The  source  of  the  low  frequency  noise  has  not  been  identified.  It  could  be  due 
to  a  slow  tilt  of  the  accelerometer  platform  brought  about  by  mechanical  creep  in  the 
suspension.  Such  a  tilt  would  modulate  the  earth’s  gravity.  Low  frequency  drifts  in  the 
oscillator  amplitude  or  the  temperature  may  aJso  be  sources  of  1//  noise.  Additionally, 
flux  creep  can  cause  drifts  in  superconducting  levitation  circuits  which  induce  1//  noise 
in  the  outputs  [9].  As  mentioned  previously,  Nb-Ti  wire  was  used  in  the  prcaent  design 
of  the  SSA  in  order  to  overcome  the  reliabilty  problems  encountered  with  the  niobium 
wire  available  at  the  time.  It  is  well  known  that  flux  creep  may  bo  significant  in  type- 
II  superconductors  such  as  Nb-Ti.  Future  work  must  concentrate  on  determining  and 
reducing  low  frequency  noise. 

Energy  Coupling  and  Fundamental  Instrument  Noise 

The  fundamental  instrument  noise  in  the  three  linear  degrees  of  freedom  in 
the  SSA  is  described  by  Cqs.  (2.229)  through  (2.232).  In  deriving  these  equations,  we 
assumed  that  all  three  linear  modes  wv  ,  •  ie  iticaJ  and  did  not  take  into  account  varia¬ 
tions  in  the  resonance  frequencies.  Wr,  tiw  .lality  factors,  Qr.  and  the  energy  coupling 
coefficients,  Replacing  (ir  by  Pry,  and  /?r,,  and  Qr  by  Qry  Qvy,  and  Qry  and 
substituting  Wr,,  i^ry,  and  in  place  of  the  generic  uir,  we  find 
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where  (ir^,  0ryi  and  0r,  are,  from  Eq.  (2.228), 

(4.75) 

(4.76) 

(4.77) 

Substituting  the  sensing  currents  along  with  other  parameiers  into  these  equa¬ 
tions,  we  find 


0r. 

=  3.3  xlO'*^, 

(4.78) 

=  2.0xl0*^ 

(4.79) 

0r. 

=  1.1  X  10~®. 

(4.80) 

These  couplings  are  extremely  small  and  explain  why  we  did  not  have  to  consider  the 
contribution  of  the  sensing  current  to  the  linear  resonance  frequencies  in  the  previous 
section.  These  small  couplings  severely  limit  the  sensitivity  of  the  accelerometer.  In  the 
new  design,  these  coupling  coefficients  should  be  larger,  and  the  SSA  should  approach 
its  potential  sensitivity. 

The  white  noise  level  normally  observed  from  the  Quantum  Design  RF  SQUIDs 
is  IQ—*  $o/V^fZt  which  is  equivalent  to  a  power  spectral  density  of  Es  = 

Substituting  this,  tlie  experimental  parameters,  and  the  values  for  the  energy  coupling 
coefficients  into  expressions  foi  the  power  spectral  density,  Eqs.  (4.72)  through  (4.74). 
we  find 
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=  1.6  X  lO-^g/yiD.  (4.82) 

Ej/2  =  5.9  X  10“’°g/v^-  (4.83) 


At  frequencies  above  approximately  40  Hz,  the  signal  due  to  seismic  noise  drops  below 
the  white  noise  level  of  the  SQUID  amplifier,  so  th's  level  ran  be  observed  directly. 
Dividing  by  the  measured  sen.sitivities,  Eqs.  (4.60)  through  (4.62),  we  obtain 


=  3.2  X  10"'‘°g/>/TD, 


(4.84) 
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(4.86) 


=  l.OxlO'Vv^- 
Pl[^  =  6.4xlO"'Vv^. 

These  axe  in  reasonable  agreement  with  the  theoretical  noise  limits. 

4,3.4  Angular  Acceleration 

The  three  angular  acceleration  circuits  are  calibrated  by  driving  the  shaker 
about  the  vertical  axis  and  comparing  the  three  angular  signals  from  the  SSA  to  the  rms 
amplitude  of  the  angular  motion. 

The  dy,  and  6^  circuits  were  driven  with  currents  of  25.5,  13.3  mA,  and  13.0 
mA,  respectively.  Multiplying  tho.se  currents  by  the  three  corresponding  gains  in  Table 
4.3,  we  find 


=  484  mA,  (4.87) 

=  146  mA,  (4.88) 

=  110  mA.  (4.89) 

The  experimental  sensitivities  for  dxi  ^y,  ^.nd  9^  at  these  drive  currents  were 

measured  to  be 


=  2.1  X  10^  4>o/rad/sec^, 

(4.90) 

=  6.0  X  10^  *o/rad/sec^. 

(4.91) 

=  4.7  X  10^  $o/rad/sec^. 

(4.92) 

The  theoretical  transfer  functions  for  the  three  angular  degrees  of  freedom  are, 
from  Eqs.  (2.198)  through  (2.200), 
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Since  actual  angular  rc,sonance  frequencies  differ  slightly,  ujg  must  be  replaced  by  ug^, 
log^,  and  ug^  in  each  of  the  three  equations.  Likewise,  Ig  must  be  replaced  by  Ig^,  Ig^. 
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and  Making  these  substitution  and  taking  the  limit  where  u  <  we  find 
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Substituting  the  values  of  sensing  currents  and  other  parameters  into  these  equations, 
we  find 


=:  2.1  X  10^ 'I>o/rad/sec^ 

(4.99) 

~  6.2  X  10^  4>o/rad/scc‘. 

(4.100) 
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These  values  are  in  good  agreement  with  the  experimental  sensitivities. 

The  low-frequency  spectra  of  seismic  noise,  as  measured  by  the  Oi,  dy,  and  9. 
sensing  circuits,  are  shown  in  Fig.  4.6.  The  figure  shows  the  fundamental  and  harmonics 
of  the  angular  calibration  signal  at  0.4  Hz.  The  four  other  peaks  are  due  to  modes  of 
the  suspension,  such  as  the  pendlum  mode  at  0.36  Hz  and  the  vertical  spring  mode  at 
1.3  Hz. 


Energy  Coupling  and  Fundamental  Instrument  Noise 

The  fundamental  instrument  noise  in  the  three  angular  degrees  of  freedom  in 
the  SSA  is  described  by  Eq.  (2.233).  In  deriving  this  equation,  we  assumed  that  all  three 
angular  modes  were  identical  and  did  not  take  into  account  variations  in  the  resonance 
frequencies,  ug,  and  the  energy  coupling  coefficients,  fig.  Replacing  3g  by  3g^,  3g^,  and 
l3g^,  and  substituting  tug^,  and  ug^  in  place  of  the  generic 
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Figure  4.6.  The  low- frequency  spectra  of  seismic  noise  as  measured  hy 
the  Bx,  By,  and  B^  .sensing  circuits. 


where  (5^^,  and  are,  from  Eq.  (2.234), 
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Substituting  the  experimental  values  of  the  sensing  currents  along  with  oilier 
parameters  into  these  equations,  we  find 


lie. 

=  2.0  X  10-^ 

(4.108) 

fiOy 

=  1.8  xlO"*^, 

(4.109) 

=  8.7  X  10~^ 

(4.110) 

These  couplings  are  extremely  small  and  explain  why  we  did  not  have  to  consider  the 
contribution  of  the  sensing  currents  to  the  angular  resonance  frequencies  in  the  previou.s 
section.  These  smaU  couplings  severely  limit  the  angular  acceleration  sensitivity  of  the 
present  model  of  the  S3A,  as  in  the  case  of  the  linear  accelerations.  The  couplings  of  the 
new  design  should  be  much  larger,  and  the  SSA  should  approach  its  potential  sensitivity. 

Substituting  Es  =  10“^*  J/Hz  for  the  Quantum  Design  RF  SQUIDs,  the  exper¬ 
imental  parameters,  and  the  values  for  the  energy  coupling  coefficients  into  expressions 
for  the  power  spectral  density,  Eqs.  (4.102)  through  (4.104),  we  find 
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e  SQUID  amplifier  can  be  observed  directly 

at  frequen- 

cies  above  40  Hz.  Dividing  this 

noise  level  by  the  iiioa.sured  semsitivities,  1 

Fqs.  (  1.90) 

through  (4.92),  we  find 
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(4.116) 

These  are  in  reasonable  agreement  with  the  theoretical  noise  limits. 


4.4  Future  Improvements 


The  present  design  of  the  SSA  works,  and  the  ineasiireinents  of  its  parameters 
agree  quite  well  wdth  theory.  Analysis  of  the  experimental  results  [)oints  out  several  areas 
in  which  improvement  cati  be  made:  1)  The  minimum  detectable  signal  can  potentiallv 
be  increased  several  orders  of  magnitude  by  increasing  the  electro-mechanical  coupling, 
3.  2)  The  sources  of  the  1//  noise  must  be  determined  and  reduced,  ‘i)  The  device  must 
be  operated  in  the  closed-loop  mode  to  increa.se  the  maximum  detectable  acceleration 
signal  and  reduce  nonlinearity. 

The  electromechanical  coupling  can  he  most  easily  increased  by  raising  the 
bridge  drive  current,  T.  Ilowevi'r,  the  misbaJance  in  the  resistive  component  of  the 
bridge  circuit  imposes  a  limit  on  the  current,  as  discussed  in  Section  •1..3.2.  In  order  to 
achieve  better  matching  between  the  coils,  it  is  necessary  to  use  single  layer  coils  for 
both  levitation  and  sensing.  However,  this  reduces  the  inductance  of  the  levitation  coils 
and  so  increases  the  current  required  for  levitation.  One  way  to  avoid  approaching  the 
critical  current  of  the  wire  is  to  use  a  single  large  roil  on  each  active  fare  of  the  roil 
form  in  place  of  the  seperate  sensing  and  levitation  coils.  This  configuration  rerpiires 
the  sensing/levitation  circuit  for  each  axis  to  have  a  separate  SC^'Hl),  but  it  has  the 
advantage  that  tiie  single  coils  are  much  easier  to  fabricate.  Furthermore,  the  single 
coil  will  have  a  much  larger  value  of  A^:  approximately  5.7  times  the  present  value.  It 
is  expected  that  with  these  improvements,  a  bridge  drive  current  in  excess  of  1  Amp 
will  be  attained.  Assuming  the  re.sonarice  frequency  is  kept  the  same  as  that  of  the 
present  model,  3  should  be  approximately  10'^  time  larger  than  current  values,  implying 
a  minimum  detectable  signal  10*  times  lower. 

As  discussed  briefly  in  Section  '1..3.3,  <  .ic  likely  source  of  1//  noise  is  flux  creep 
in  Nb-Ti.  F'lux  creep  should  be  greatly  reduced  by  using  Nb  in  all  the  superconducting 
circuits.  Switching  back  to  Nb  at  this  point  is  reasonable  because  of  a  better  understand¬ 
ing  of  how  the  processing  of  the  Nb  wire  affects  its  mechanical  properties.  In  addition, 
using  only  a  single  sot  of  coils  for  both  levitation  and  sensing  siitiplifics  the  overall  design. 

Like  most  conventional  accelerometers,  the  SSA  should  operate  best  in  a  feed¬ 
back  loop:  The  six  demodulated  outputs  are  fed  into  a  controller  which  applies  current 
to  the  six  levitation  circuits  to  maintain  the  balance  of  the  bridges.  In  this  way,  the 
effective  bandwidth  of  the  device  can  be  pushed  to  several  times  the  natural  resonance 


100 


frequency.  The  controller  also  provides  damping.  Because  the  o[)ei:  loop  system  is  lightly 
dii.mped,  a  large  acceleration  impulse  or  strong  acceleration  noise  oflcn  causes  a  signal 
at  the  resonance  frequency  large  enough  to  excec  1  the  slew  rate  I'f  the  SQl'ID  am{)li 
fier.  By  providing  damping,  the  maximum  readable  signal  is  incrca.sed.  In  the  [)resf'nt 
circuit,  all  the  sensing  bridges  are  directly  connected,  so  the  spacing  between  the  rarner 
frequencies  limits  the  bandwidth  allotted  to  each  axis,  and  thus  the  iltimato  efferiive 
bandwidth  attainable  using  a  controller.  By  using  separate  SQTlDs  for  each  axis,  t;  <' 
ultimate  bandwidth  of  the  accelerometer  should  be  limited  by  the  carrier  frequency  if 
coupling  between  the  circuits  through  stray  inductatices  can  be  made  negligible. 

A  new  version  of  the  accelerometer  has  boon  designed  incc)rporating  the  above 
changes.  The  Model  II  SS.\  will  have  only  21  coils  connected  in  six  independent  bridge 
circuits  that  will  provide  both  levitation  and  sensing.  Each  bridge  will  have  its  own 
SQUID  amplifier.  A  simple  BID  (r-oportional,  integral.  diff(>rent iai )  controller  h;u^  re¬ 
cently  been  built  and  te.stod  on  the  present  version  of  the  SS.\.  Ucing  the  results  of  thes<' 
tests,  an  improved  controller  will  be  built  for  the  .Model  11  SS.\.  With  these  imj)rove- 
ments,  substantial  improveitient  in  sensitivity,  especially  in  the  crucial  low  fft'ciumu >' 
region,  should  be  attained. 
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